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PREFACE TO THIRD EDITION 


In preparing this third edition it has not been considered 
necessary to make any change in the lay-out of the text, 
the general form of which seems to have served the needs 
of the physicists and engineers for whom the book was 
written. 

Opportunity has been taken, however* to add some 
further notes on vector notation (p. 117), on the more 
formal aspect of the del operator (p. 113) and on Maxwell's 
equations (p, 118). At the same time a number of minor 
alterations and corrections have been made, I am indebted 
to Professor S. Chapman* F.R,S,, and to Professor G* 
Cook* F,R, 3 *, for kindly criticisms and helpful suggestions 
which have enabled me to make the text clearer and more 
accurate. For certain historical facts I am indebted to 
the late Dr, John MeWhan* 

Glasgow 
Jan, 1945. 


PREFACE 

Vector Analysis is the natural means of expression for the 
three-dimensionai problems of physics and engineering, 
because its conciseness and freedom from mathematical 
detail enable the relationships between the various physical 
quantities to be kept clearly in view. Since the pioneer 
work of Gibbs and of Heaviside an increasing number of 
text-books and scientific papers on physical and technical 
subjects have made use of vector methods, until it has now 
become almost essential for any advanced worker in these 
sciences to have some knowledge of vector analysis. Much 
good work can be done with the aid of a very few elementary 
principles. It is the object of this monograph to give an 
introduction to these principles and to explain them from a 
physical standpoint, so that they may be easily available to 
the busy physicist or engineer approaching the subject for 
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the first time. Such workers are usually so much occupied 
by their major task as to lack the time necessary to enable 
them to seek out such principles as they need to use from 
the more comprehensive treatises which aim at mathe¬ 
matical completeness. The monograph is not, therefore, 
intended for the reader with purely mathematical interests, 
whose more rigorous and systematic requirements are fully 
satisfied elsewhere. 

For these reasons the outlook adopted is almost entirely 
physical; geometrical matters and questions of an exclu¬ 
sively mathematical interest are limited to essentials. 
Formal proofs of invariance and conditions of continuity 
in vector processes are replaced by an appeal to physical 
intuition. Purely analytical topics of an advanced kind, 
such as Green’s theorem, are justifiably omitted from such 
an introductory treatment as this monograph aims to pre¬ 
sent. The reader who later wishes to amend his knowledge 
in these and other respects will find ample material in the 
standard books listed in the Bibliography, 

The reader’s attention is drawn in particular to two 
features. First, the use of trimetric projection by Gough's 
method {Engineerings voh 143, p. 458, 1937) for certain of 
the three-dimensional diagrams, a purpose for which this 
method is admirably suited. Second, the inclusion of a 
chapter giving a brief sketch of the elementary properties 
of tensors and dyadics in their relation to vectors. This is 
a subject which usually puzzles and often repels physics 
and engineering students because of the abstract mathe¬ 
matical way in which it is generally brought to their notice. 

The monograph is based on a course of lectures given 
a few years ago to post-graduate electrical engineering 
students in the Polytechnic Institute of Brooklyn, New 
York. I am particularly grateful to my colleague Dr. A, J. 
Small for his valued assistance in reading the manuscript 
and proofs, and for numerous suggestions. 

Glasgow 
June 1938 
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CHAPTER I 


DEFINITIONS, ELEMENTS OF VECTOR 
ALGEBRA 

I. Scalar and Vector Quantities. Physical quanti¬ 
ties arc divided into two main classes* each with charac¬ 
teristic properties and an appropriate algebra. 

Scalar quantities have magnitude only and do not involve 
direction. Typical scalar quantities are mass, volume, 
density, temperature, electric potential, charge, &c. The 
complete specification of a scalar quantity requires (i) a 
unit of the same kind and (ii) a number stating how many 
times the unit is contained in the quantity. For example, 
to express the mass of a given body we require to know 
whether the unit is the pound, the ton, the gramme, 
and also how many of the chosen units represent the given 
mass. Scalar quantities are manipulated by applying the 
rules of ordinary algebra to their numerical magnitudes ; 
for this reason the algebra of ordinary positive and negative 
numbers is often called scalar algebra and the numbers 
themselves scalars. 

Vecior quart titles have m agn i tu de and d irecti on, Fam i 1 iar 
examples are displacement, force, velocity, acceleration, 
stress, electric force, magnetic induction, &c, A vector 
quantity requires for its specification (i) a unit of the same 
kind, disregarding direction, (ii) a number giving the mag¬ 
nitude of the quantity in terms of this unit and (iii) a 
statement of direction. For example, the velocity of a 
moving body is stated by saying (i) that the unit is miles 
per hour, kilometres per second, &c, ; (ii) how many of 
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the chosen units express the magnitude of the velocity 
and (iri) the sense in which the velocity is directed, e.g. 
due north. The combination of conditions (ii) and (iii) 
constitutes the geometric conception of a directed magni¬ 
tude or vector, quite independent of the kind of vector 
quantity specified by the unit. The directional element 
will prevent the manipulation of vectors by the simple 
numerical algebra applying to scalars ■ it is essential, there¬ 
fore, to devise a vector algebra by means of which vectors 
may be handled in a way consistent with the physical 
problems in which vector quantities occur. The laws of 
vector algebra differ in several important respects from 
those of scalar algebra, as will be seen later. 

Vector quantities, though the simplest, are not the only non- 
scalar quantities found in physics and engineering. In the theory 
of inhomogeneous strain, in the transformations of space involved 
in the theory of relativity and in the general theory of electro¬ 
magnetic machinery, to mention only a few examples, more com¬ 
plex non-scalars appear. These are the linear vector function, 
the dyadic and the tensor. Each has its own algebra, differing 
radically from the scalar laws and more general in form than 
vector algebra; indeed, the particular quantity is characterized 
by its algebra. To manipulate these quantities with conciseness 
and ease, higher algebraical methods, multiple algebras and the 
theory of matrices are widely used. In this elementary mono- 
grt^ph it is not possible to do more than sketch the simplest 
properties of these higher non-scalars and their relation to ordinary 
vectors ; the reader will find further developments in the special 
treatises referred to in the Bibliography. 

Scalars arc manipulated by six elementary operations, 
namely, addition and subtraction, multiplication and 
division, involution and evolution. In each pair the 
second operation is the inverse of the first. Also, multi¬ 
plication is extended addition, and involution (i.e. raising 
a number to a power) is extended multiplication ; like¬ 
wise, subtraction, division and evolution are successively 
extended inverse operations. Ail these operations are 
performed by means of certain laws, which will now be 
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stated for addition and multiplication ; they are readily 
extended to the other operations and the reader will find 
formal proofs in any good text-book of algebra : 

(i) The law of association can be written as 

a + (ft + c) - (a -f ft) + c t 
and a x (ft x c) = (a * ft) x c; 

i.e. the result of adding a to the sum of ft and c is the same 
as adding the sum of a and ft to c ; with a similar state¬ 
ment for products. In other words, brackets are unneces¬ 
sary in continuous sums and products. 

(ii) The laic of commutation which states that additions 
and multiplications can be made in any order; i.e. 

a + ft - ft + a f 
and a x ft *- ft x a, 

(iii) The law of distribution, stating that in compound 
expressions involving multiplication and addition the result 
is equivalent to the sum of individual products taken term 
by term, as in the rules 

a x (ft + c) - a x ft + a x c w 
and (fl + ft) xc-tfxc + ftxf. 

The three laws are extended to the inverse operations 
by the rule of signs, and additional rules for the manipula¬ 
tion of indices in involution are also provided. Finally, 
a fourth law, fundamental in the logical discussion of inverse 
operations, is now usually stated, namely, if a x ft ~ o, 
then either a - o or ft » o. 

While these scalar operations and their laws are familiar 
to all as the working processes of arithmetic and algebra, 
so familiar indeed as to be regarded as self-evident, it is 
quite essential to state them in this formal way since they 
characterize scalar algebra. The essential feature of the 
algebras of vectors, tensors, matrices, &c., is that they all 
violate in some way the commutative law for multiplica¬ 
tion and do not necessarily satisfy the condition that 
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either factor in a product is zero if the result is zero. The 
physical meaning of these differences between scalar and 
vector algebra will be explained later* 

2 . Graphical Representation of Vectors* Since a 
vector is the result of abstracting magnitude and direction 
from a vector quantity, independently of the nature of the 
quantity concerned, we may represent the vector graphi¬ 
cally by a line CM pointing in the direction from O to A f 
as in Fig* i. The magnitude of the vector is given to a 



Fig. i .-“Definition of Equal Vectors 


convenient scale by the length of the line ; the direction 
in space is indicated by an arrow-head marked on the line* 
In the algebra of such vectors it is the custom to distinguish 
them by the use of a distinctive symbol* Thus, the vector 
0*4 is denoted by the Clarendon letter V ; its scalar magni¬ 
tude is stated by the corresponding italic letter F. A 
vector is unchanged by pure displacement; it follows, 
therefore, that all parallel vectors of the same length and 
direction are equal ; a vector is not localized. Thus in 
Fig* i the vectors V, V lf V*, V> are all equal, a fact stated 
by the relation 

v-v.-v.-v* 

in which the meaning of the sign « is extended to include 
equality of size and similarity of direction. Reversing the 
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arrow-head changes the sign of a vector, such as V 4 in 
Fig* i; then 

V, - - V 

indicates equality of size and opposition of direction. 

This definition restricts the meaning of a vector to the repre¬ 
sentation of size and direction only. Consequently, a vector 
cannot represent completely any vector quantity the effect of which 
is changed by simple displacement; many quantities are of this 
kind, most notably force. It is clearly insufficient to specify a 
force by its magnitude and direction alone, i.e* by its vector. 
If the force acts on a rigid body the line of action must also be 
stated, since a change in the line of action alters the torque acting 
on the body. In the case of a deformable body the point of 
action of a force must be specified. These examples are given 
as a warning to the reader not to be tempted to attribute to a 
vector any properties, other than those of size and direction, which 
properly belong to the quantity that the vector is representing to 
scale merely in these particulars* In other words, vector analysis 
must always be used with an eye upon the physical conditions of 
the problem to which it is applied* 

It is geometrically obvious that multiplication of a 
vector by a scalar factor S results in a vector 5 times the 
size of the first and in the same direction* in particular, 
if v be a vector of unit size in the direction of V, 

V ~ Fv.(r*i) 

Distinction is often made, particularly in more advanced 
work, between two essentially different classes of quantities 
represented by vectors* The first class includes quantities 
such as force, displacement, velocity, &c*, in which the 
vector is drawn in the direction of the quantity concerned ; 
mere linear action in a particular direction is involved and 
the vector is called a polar vector. The second class 
includes quantities such as angular velocity, angular 
acceleration, See*, in which rotary action of some kind takes 
place about an axis* Here the vector is drawn parallel 
to the axis about which the quantity acts ; the length of 
the vector gives the magnitude of the quantity* The 
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direction of the vector is conventionally fixed by the rule 
of the right-handed screw, i.e. when sighting along the 
direction shown by the arrow-head ro ™ !* 
clockwise. Such vectors are named 
distinction is a physical one, there being little dil 
ence in the mathematical treatment of the two classes o 
quantities. Indeed, it is difficult in many cases, particu¬ 
larly in dynamics and electromagnetic theory, to assig 
vector quantities with certainty to classes in this way 
3 . Addition and Subtraction of Vectors. Consider 
two'vectors A and B, shown in Fig. 2(a), representing or 



(«) 




F»- *-« Vector V ' C '" 

example two successive displacements of a point. Their 
ioffiTeffect, sum or resultant is obtained by setting off the 
vector B at the end of A and drawing the vector V joining 
Ihe beginning of A to the end of B. The same result rs 
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obtained by starting with B and attaching A to it m a 
similar way, as shown by the dotted lines* Then 

V = A + B»B + A, , , , . (i.2) 

the + sign being understood to mean addition in this 
geometric sense* Thus the sum of two vectors U the 
diagonal of the parallelogram of which the vectors are the 
sides ; such a sum is commutative, i.e. independent of 
which vector is taken first. Subtraction of one vector from 
another immediately follows ; the vector B to be sub¬ 
tracted is reversed and the rule for addition applied, as 
illustrated by Fig* 2 (b) r in which V- A + (— B) - A — B. 

A useful geometrical application is shown in Fig* 2(c), 
where A and B are given vectors. If kk a variable para¬ 
meter ranging between + co and - oo, the sum of A and 
kB is the radius vector from O, namely, 

r -A + AB. 

As k varies, the locus of P will be a straight line through A 
parallel to B, and the expression is the vector form of the 
equation for a straight line. 

The reader must not suppose that all quantities that may be 
represented by a directed line are necessarily vector quantities ; 
the crucial test is whether they follow the parallelogram law of 
addition or not* For example, finite rotation of a rigid body 
about an axis can be represented by a line, of a length propor¬ 
tional to the angle of rotation, drawn in the direction of the axis, 
A second finite rotation about another axis inclined to the first may 
be similarly represented. But their resultant effect cannot be 
found by compounding the two finite rotations by the law for the 
vector sum, as the reader may easily verify in a simple case such 
sa the motion of a hall on a plane. Finite rotations are not, 
therefore, capable of representation by vectors ; they require the 
use of tensors and have much more complex laws of manipulation. 
Infinitesimal rotations, angular velocities and accelerations are, 
on the other hand, vector quantities (see p, 46}. 

With more than two vectors, such as A, B, C in Fig* 3, 
the sum is obtained by laying off the vectors successively 
2 















« 


VECTOR ANALYSIS 


end to end ; then the sum is the dosing side V of the 
polygon of which the other sides are A, B and C + But it 
is geometrically obvious that we could start by adding any 



F (G . ^—Law of Association for Vector Addition 

pair of vectors by the parallelogram rule, and to their 
resultant add the third vector. Hence 
V-A+B+C«(A+B)+C-A+(B+C)-(A+C)+B. . {1*3) 

The sum of any number of vectors is, therefore, associative, 
j.e. the vectors may be added in any desired manner. The 
vectors are not necessarily all in the same plane ; this is 
more clearly shown by Fig. 5* 

4* Components of a Vector* By reversing the 
process of addition it is clear that any vector V can be 
decomposed into the sum of n vectors, of which n - 1 are 
arbitrary and the last one closes the polygon* In general 
the vectors are not eoplanar and the polygon is not a 
plane figure* 

The most useful instance is the decomposition of a 
vector into component vectors along the three orthogonal 
axes of cartesian co-ordinates. Fig. 4 {d) shows the usual 
right-handed system, in which the relative positive direc¬ 
tions along the axes are chosen such that if O'X is turned 
toward O' Y about O'Z through the smaller angle, a right- 
handed screw would advance along the positive direction 
of O'Z ; similar relations hold for the axes O' Y and 0 *X t 
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always keeping X $ Y and Z in cyclic order. Let the origin 
O f be one extremity of a vector V ; draw a rectangular 
parallelepiped with the three edges which meet at O' 


r 




Fig, 4,—Cartesian Components of a Vector 


lying along the axes and such that V is the diagonal from 
O' through the solid figure. If V„ \ y , V. are the 
vector intercepts along the axes of X, Y and Z respectively, 
then 


v - V, + V„ + V,. . 


• • (1+) 
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Now let i, J t k denote vectors of unit magnitude along 
the axes of X , Y and Z respectively ; these are the unit 
vectors along the axes and are frequently used in later work. 
Then if V xt V yl V f are the sizes of the vectors V*, V*, V„ 

V, - VJ, V v - VJ and V f - 

from Equation (l*l) ; whence 

V - VJL + VJ + V 2 k .(i* 5 ) 

In discussing the properties of a vector field we are 
concerned with the magnitude and direction of a vector 
quantity at ever) 1 point in space. In rectangular co¬ 
ordinates with origin O let a point in space, such as 0 in 
Fig* 4(&), be specified by the radius vector 

r ~ xl + yi + ~k* 

Since a vector is not affected by parallel displacement, 
Equation (1.5) is still the relation between the components 
of a vector and the vector itself. Note, however, that as 
we are dealing with a vector field, V is a function of posi¬ 
tion ; that is, V and its components V„ V yi V t are func¬ 
tions of the variables y t z specifying the position of O 
relative to 0* 

Vector quantities and their vectors are physically inde¬ 
pendent of any system of co-ordinates by means of which 
they may be expressed ; they are said to be invariant. 
Co-ordinate axes are unnecessary, therefore, in vector 
analysis. The resolution of a vector into components, 
given by Equation 1.5, is the connecting link between 
vector notation and that of ordinary co-ordinate geometry 
with which the reader is already familiar. Since a vector has 
three axial components, any vector equation relating vectors 
is the equivalent of three scalar cartesian equations relating 
their components* This conciseness is the great advantage 
of vector methods in three-dimensional problems* The 
fact that vector quantities are expressed directly by a nota¬ 
tion natural to themselves, instead of by a long and artificial 
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system of scalar relationships, enables their essential 
physical nature to be kept dearly in view* The resolution 
into cartesian components is often useful in the proof of 
theorems used in vector analysis ; in the words of Heavi¬ 
side, 4 When in doubt and difficulty, fly to i, j, k \ advice 
we shall frequently follow. 

As a simple example find the sum of vectors A, B, C. In 
terms of their components, 

A - AJ + A v j + AX B ~ BJ + BJ + B £ k, 

C - C x i + CJ + C 2 k. 

The components along the X axis add directly ; like¬ 
wise those along the Y and Z axes* Then, as shown in 
Fig- 5* 


Y 



Fig* 5*—Cartesian Components of Vector Addition 

V * VJ. + VJ + F*k - A + B + C 
- (A* + B x + C £ )i + (A, + By + C y )i 

+ (A; + B t + CJ k . * (1*6) 

Since the scalar magnitudes in the brackets obey the 
ordinary laws of algebra, it follows that vector sums (and 
differences) obey these laws also. The component of the 
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resultant in any axial direction is the sum of the individual 
components in that direction ; as is geometrically obvious 
from Fig. 5, 

V x - A % + Vy"*Ay + B v + C v and 

V t - A 2 + B z + C a . 

5 , Scalar and Vector Fields, A physical quantity 
can be expressed as a continuous function of the position 
of a point in a region of space ; such a function is called a 
point-function and the region in which it specifies the 
physical quantity is known as a field . Fields are of two 
main kinds, scalar and vector, according to the nature of 
the quantity concerned. 

A typical scalar field , such as the distribution of tem¬ 
perature, density, electric potential or of any other non- 
directed quantity, is represented by a continuous scalar 
function giving the value of the quantity at each point. 
Such a function does not undergo any abrupt change of 
magnitude in passing from any point to another close to 
it, a condition satisfied in all practical cases. The held can 
he mapped graphically by a series of surfaces—such as iso¬ 
thermal, equi-density or equipotential surfaces—upon each 
of which the scalar has a definite constant value. Such 
surfaces, called equal or level surfaces, are conveniently 
chosen so that in passing from one to the next a constant 
arbitrary difference is made between the scalars which 
characterize them. It is evident that the level surfaces 
must lie one within the other and cannot cut ; for if two 
such surfaces could intersect, the scalar values corre¬ 
sponding to both must hold along their common line, 
which is contrary to our definition. Hence scalar point - 
functions are single valued or uniform at every point. 

A typical vector field , such as the distribution of velocity 
in a fluid or of electric or magnetic field strength, is repre¬ 
sented at every point by a continuous vector function. At 
any given point the function is specified by a vector ot 
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definite magnitude and direction, both of which change 
continuously from point to point throughout the field 
region. Starting at any arbitrary place, proceed an in¬ 
finitesimal distance in the direction of the vector at that 
place, arriving at a closely-neighbouring point. Proceed¬ 
ing thence in a similar way, we shall trace out a curved 
line, the tangent to which at any point gives the direction 
of the vector thereat; such a curve is a vector line, line of 
flow or flux line. To represent the magnitude of the vector, 
at any point on a flux line draw a very small surface per¬ 
pendicular thereto and choose a number of points per unit 
area upon this surface numerically equal to the magnitude 
of the vector. Through each of these points flux lines 
can be drawn. The field is then mapped out by flux lines. 
The direction of the lines is that of the vector function ; 
their density, represented by the number of them crossing 
per unit area perpendicular to their direction, is a measure 
of the magnitude of the vector. It is dear that lines of 
flow cannot intersect, since this would involve indefinite 
direction of the vector at the point where they cut ; vector 
point-functions must also be single valued at every point. 

The physical properties of scalar and vector fields will 
be considered in greater detail in Chapters VII and VIII 
respectively. 












CHAPTER II 

THE PRODUCTS OF VECTORS 


1. General. The ordinary idea of a product in scalar 
algebra* the mere multiplication of a scalar magnitude, 
cannot apply to vectors because of their directional proper¬ 
ties ; nor is it possible to decide by deductive reasoning 
what form the product of two vectors should take. Since 
vectors have their origin m physical problems* definitions 
for the products of vectors must be devised that will be 
consistent with the way in which such products occur in 
applications to physical science. 

As a simple example consider vectors which represent 
a force F and a linear displacement d, their magnitudes 
being F and d and the angle between their directions 0 . 
Products of these two quantities occur in two ways. First, 
the work done by the force is Fdcos0, a scalar quantity, 
known as the scalar product of F and d, The couple 
exerted by the force has a magnitude Fd sin 0 and acts 
about an axis perpendicular to the plane containing F 
and d, It is shown by a vector, called the vector product 
of F and d, drawn along the axis of the couple in a con¬ 
ventional sense to be defined later* We arc led, therefore, 
to define two sorts of products, namely* the scalar product 
and the vector product, 

2 . The Scalar Product- The scalar product of two 
vectors A and B is defined as the product of the magnitudes 
of the vectors and the cosine of the angle between their 
directions. As shown by Fig. 6, this is the same as the 
product of the size of one vector with the component of 
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the other in the direction of the first* and this process is 
clearly commutative* since cos 0 cos {- 0), i.e. indepen¬ 
dent of the order of the factors. The scalar product will 



Fig. 6 .—Scalar Product of Two Vector* 


be denoted by interposing a dot between the vectors ; 
then 

A B - B*A AB cos 0 . , . (2.1) 

Thus* by definition* the scalar product is commutative. 
On p. 3 it has been pointed out that the product of two 
scalars vanishes when either is zero. The scalar product 
of vectors also satisfies this condition* but will vanish in 
addition when the vectors are at right angles. This is the 
first particular in which the laws of vector algebra differ 
from those of scalar algebra. When two vectors are 
perpendicular* therefore* 

A'B “ o* • * • # • (2.2) 

and when they are parallel, 

A B - AB .. (2,3) 

If B - A t the scalar product of a vector with itself, called 
its self-product, is 

A-A-A . . . .( 2 . 4 ) 

hence to find the size of a vector calculate the square root 
of its self product. 











VECTOR ANALYSIS 


16 

Applying these properties to the unit vectors I, J, k, 
which are mutually perpendicular, 

ij = jk - k i - o, .... {2.5) 
and . . . .(2.6) 

Now examine the scalar product of a vector A with the 
sum of two others, B and C, illustrated in Fig. y(a ); the 



(a) 



<*) 


Fig* 7* — (a) Law of Distribution for Scalar Products ; (6) Vector 
Equation of a Plane 

three vectors are not necessarily in the same plane. Since 
the resolute of B + C on A is equal to the sum of the 
resolutes of B and C, i,e. OC - OB + BC , 

A-(B + C) - A(OC) « A(OB + BC) - A*B + AC. (2*7) 

Hence the scalar product is distributive for addition, as 
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is geometrically obvious* By an extension of this process 
it is easy to show that 
(A+B-k , .)<N+0+..,) - 

A-N+A* 0 +. t( B-N+B-0+... . . ( 3 . 8 ) 

A particular case is of geometrical interest. In Fig, j(h) 
let p be the perpendicular from O upon a plane which 
passes through the extremity of a given vector c. If r 
is the radius vector to any other point in the plane, then 
it is obvious that the resolutes of c and r upon p must 
be equal, he. 

P -r - p c or p (r - c) - o 

is the equation of a plane perpendicular to p through the 
end of c. From Equation 2,2, {r - c) must lie in the plane 
normally to p, as is obvious. 

Referring again to Fig. 2(a), the scalar product enables 
the usual formula for the resultant of two vectors to be 
readily found. Their sum is 

V - A + B* 

Take the self product on both sides ; then from Equations 
2.4 and 2.8, 

V*V - V 2 - (A + B)*(A + B)« A A + 2A B + B B, 
or V 2 - A % *f 2 AB cos 0 + 

The rectangular components of a vector are found at 
once by taking its scalar product with the unit vectors, as 
is evident from Fig, 4 ; then 

V i - V xt V j - V y and V k - V t * . . (2-9) 

The scalar product assumes an important form in 
rectangular co-ordinates. Writing from Equation 1.5, 

A - AJ + A J + A ,k and B - B x i + B y j + B, k. 
Equations 2.5, 2.6 and 2.8 give 
A-B - (AJi + A y j + A*k){BJ + By j + B,k) 

“ A X B X + A y By + A z B $m # (2.10) 
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That is, the scalar product of two vectors is the sum of 
the products of their components along each of the co¬ 
ordinate axes. If, for example, A is a force and B is a 
displacement, the total work done is equal to the sum of the 
works done by the components of force and the correspond¬ 
ing components of displacement, as is physically obvious. 

3 . Line and Surface Integrals as Scalar Products. 
In Fig. 8, let OA be any curve drawn in a vector field and 



dl an element of length along it at any point P. Let V 
denote the vector at Pin a direction making an angle 0 with 
that of the length clement. Then is the product of 
the length of the clement and the component of V in its 
direction, he. tangentially to the curve ; thus 
Vdl - Fcos 0 dl 

If V varies in magnitude and direction from point to point 
along the curve, the integral 

pV-dl- rVcoafldZ * - (2.xi) 

Jo Jo 

is defined as the line integral of V along the curve OA . 

Such integrals are of very frequent occurrence. For 
example, if V is a force and dl an element of the path of a 
particle along any curve, the line integral denotes the work 
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done in displacing the particle from O to A. Again, if V 
is the electric field strength, i,e, the force on unit charge in 
an electric field, then the line integral expresses the 
potential difference between the points O and A, Further, 
if V is the velocity at any point in a fluid and the integral 
be taken round a closed curve, the integral is called the 
circulation of the fluid. In an electromagnetic field the 
line integral of the electric force round a closed path is the 
electromotive force in the path ; in a purely electrostatic 
field the ejn.f, is zero. 

Consider now any element of area ds upon a surface 
drawn in a vector field, and let V be the value of the vector 
at the middle of the element. Draw the positive* normal 
of unit length n on the element and let 0 be the angle 
between n and V, as in Fig. 9* Then the component of 



Fig. g.—Normal Surface Integral of a Vector 

V perpendicular to ds is V n « V cos 0 and the flux of V 
through the element is V-n ds. The integral of this taken 

* If the surface is dosed, n is drawn outward from the enclosed 
volume. If the surface is unclosed, then n is drawn always on the 
aume side of the surface, the normal on the opposite side being 
negative. 
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over the surface is called the total flux or surface integral of 
V through the whole surface r, i.e. 

Il. v ' ni -1L V cos 0 ds. , . (2*12] 

The physical meaning of such surface integrals is very 
simple. Suppose V to denote the vector velocity of a 
moving fluid in which a fixed surface s is drawn. At any 
point V-n ds denotes the amount of fluid passing normally 
through the surface element ds in unit time* and this com¬ 
ponent alone must be considered since the tangential com¬ 
ponent of V necessarily contributes nothing to the flow 
through the surface element. If the flow through an ele¬ 
ment is in the sense of the positive or outward normal 
it is counted as positive ; if in the opposite sense, as on 
the element ds u It is negative. The integral of these 
normal contributions expresses the aggregate flow of fluid 
through the whole surface in unit time. If the surface is 
dosed, positive total flux diverges from the enclosed volume 
while negative flux converges upon it* Should the amount 
of flux entering by some elements be balanced by that 
leaving by other elements in such a way that the total flux 
is zero, then either there are no sources or sinks of fluid 
within the enclosed volume or their sum is zero. Similar 
ideas apply to other fluxes, e.g. of electric or magnetic 
induction, of heat, 3 cc, 

4 * The Vector Product, The vector product of two 
vectors A and B is defined as a vector having a magnitude 
equal to the product of the magnitudes of the factors and 
the sine of the angle between their directions, and a direc¬ 
tion perpendicular to the plane containing A and B, The 
sense of the product along this perpendicular is defined by 
the right-hand screw rule ; if the vector A is turned to¬ 
wards B through the smaller angle, the necessary rotation 
must be clockwise when sighting along the positive normal 
unit vector n in Fig, 10* That is, the rotation needed to 
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move A to the position of B and the positive direction of n 
are related in the same way as the rotation and translation 
of a right-handed screw ; this direction of n is that 
of the vector product. It follows from this definition that 
a change in the order of the factors in a vector product 
reverses the sign of the product* since sin f- 0) = - sin 0, 

A*B 



* 

B*A —-A*B 

Fig, io.—V ector Product of Two Vector* 

Hence the vector product of two vectors is not commutative 
and the order of the terms must be strictly maintained, 
a further difference between vector and scalar algebras. 
The vector product will be denoted by a cross between the 
factors ; then 

AxB**-BxA« AB sin 6 n, . * (2.13) 

where n is a positive norma! of unit length drawn from the 
plane containing A and B. The product vanishes not 
only when either factor is zero but also when the com¬ 
ponent vectors are parallel, again differing from scalar 
algebra* When two vectors are parallel, 

A * B - o, .... (2.14) 

and when they are perpendicular, 

A x B = AB n, * , 


* * ( 2 -* 5 ) 
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in which case the two vectors and their product are mutually 
at right angles. Applying these ideas in particular to the 
cartesian unit vectors i, j, k, 

lxi = jxj-kxk = o, . • * {2.16) 

and 

Jxj—jxl=k, j*k—kx]*-i, kxl—ixk-j. . (2.17) 

The strict cyclic order in these important results should be 

observed. . , , 

5 . Vector Area. The magnitude of a vector product 

is AB sin 0 , i.e. the area of the parallelogram with sides 
A and B and included angle 0 . Since the shape of the 
figure is not specified by the definition of the product, any 
plane area of amount AB sin 0 with its positive unit normal 
n can be taken to represent a vector product. By an 
extension of this idea we are led to the notion of vector 
area. A plane area, such as s in Fig. 11(«), can be regarded 
as possessing both magnitude and direction. Its magnitude 
is the amount of the area and its direction is that ol the 
normal to its plane. The sign to be attributed to a vector 
area is defined with reference to the order in which it is 
traced out as viewed from an external point such as O. 
If the direction of tracing is counterclockwise as seen from 
this viewpoint, the positive direction of the vector area is 
along the unit normal n, the directions of tracing and of n 
being related by the right-hand screw rule. Then s = m 
is the definition of a vector area S. Vector areas can be 
resolved or added just like other vector quantities. 

In Fig. 11(A) a tetrahedron is shown with vectors s u 
s„ Sj, s, drawn to represent the vector areas of its four 
triangular faces, the outward normal being regarded as 
positive. Resolve these areas upon any plane, e.g. one of 
the faces of the solid. Then some of the projections will 
be positive and others negative, the sum of them all being 
zero. Hence the total vector area of a tetrahedron is 
zero. Thi3 result follows at once from physical considera- 


THE PRODUCTS OF VECTORS 


23 


tions, if the tetrahedron is supposed to be drawn within a 
fluid which is in equilibrium under hydrostatic pressure. 
Each face experiences a force normal to its plane and 
proportional to its area. Since the fluid within the solid 
figure is in equilibrium with that outside it, the resultant 
of the forces on its faces is zero ; hence also is the sum of 
the vector areas, since the pressure is the same on all faces. 
r i his hydrostatic demonstration also applies to any form 
of solid figure and thus generalizes the theorem. Geo¬ 
metrically, any polyhedral figure may be divided up into 



(«) (b) 

Fio. 11. —(a) Vector Area ; (A) Vector Area of a Closed Surface 


tetrahedra ; every surface introduced into the interior of 
the polyhedron appears twice, once with a positive and 
once with a negative normal. Hence for any polyhedral 
surface the total vector area is zero. By making the faces 
vanishingly small and indefinitely increasing their number 
we approach a closed curved surface over which 

6. Application to Vector Products. Consider now 
the vector product of a vector A with the sum of two 
others, B and C. In Fig. 12 draw a triangular prism with 
its parallel edges in the direction of A and its end faces as 
triangles with sides B, C and B + C. The vector areas 
ol the triangular end faces are |B x C and JC x B, which 
3 
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cancel; the remaining vector areas are B x A, C x A 

and A x (B + C) in the sense of the outward normals, 

their sum being zero. Thus 

Ax(B + C) + BxA + CxA-0 

i.e. Ax(B+C)—BxA-CxA-AxB+AxC . . (2.18) 



Fig, i a,—Law of Distribution for Vector Product 


Hence the vector product is distributive, but the order of 
the factors must be strictly observed. By repeating the 
process it is easy to show that 

(A+B-k . .)x(N+ 0 +. . .)- 

AxN+AxO+. .. BxN+BxO+.. {2.19) 

for any number of vectors. 

These rules enable the vector product of two vectors to 
be expressed in rectangular components. Using the nota- 

i tion of Equation (1,5), 

A xB= (AJ+AJ+A.k) x (S T i +B,,i + DM) 
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on applying the rules for vector products of unit vectors 
and the distributive law, Equations 2.16, 2.17 and 2.19. 
1 his result is more easily remembered by writing it in the 
form of a determinant; thus, 


A x B - 


i 

A x 

B x 


i 


A 

A 


k 

A, 


. . . (2.20 bis) 


7 . Products of Three Vectors. The vector product 
of two vectors B and C being a vector, can give both a 
scalar and a vector product with a third vector A. There 
are, therefore, two triple products, namely A(B x C) and 
A x (B x C), which occur frequently in physical applica* 
tions. . 

i he scalar triple product has a simple interpretation 
illustrated by I'ig. 13. (B x C) is a vector normal to the 


B*C 



plane of B and C, with a magnitude equal to the area of 
the shaded parallelogram. The scalar product of A with 
(B x C) is the product of this vector area and the projec¬ 
tion of A along (BxC); i.e., A'(B x C) is the volume of 
the parallelepiped which has A, B, C for its edges. Any 
face of this solid can be taken as the base ; hence three 
equivalent expressions for the volume are 

A-(B x C) - B (C x A) - C-(A x B), . (2.21) 
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cyclic order of the factors being maintained to retain the 
volume with positive sign,* Since the order of terms 
in a scalar product is immaterial, Equation 2.1, these rela¬ 
tions are equivalent to 

(A x B)-C - {R x C) A - (C x A) B ; . (2.21) 
so that the dot and cross may be interchanged at will. 


r 



Three vectors, therefore, have six identical scalar triple 
products which may be written concisely as 


11 j k 


A x Ay A z 

B 9 By B t 


B z By Bj 

C X Cy C Z 


G z Gy C z 


. (2,21 £&) 

When three vectors lie in a plane, the volume of the 


* Scalar quantities of this kind, in which the sign depends on 
the proper cyclic arrangement of the component vector factors, 
are called by mathematicians * pscudoscalars % to distinguish them 
from true scalars, which do not change sign when the reference 
axes are changed from a right-handed to a left-handed system. 
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parallelepiped is zero ; hence the condition for vectors to 
be coplanar is that their scalar triple product vanishes. 
The vector triple product A x (B x C) can also be illus¬ 
trated geometrically, as in Fig. 14, The vector (B x C) 
is normal to the plane containing B and C, Likewise, the 
vector A x (B x C) is normal to the plane containing 
A and (B x C), he. in the same plane as B and C« Take 
the X axis along B, the Y axis at right angles thereto in the 
plane of B and C, and the Z axis along (B x C). Then 
in terms of their axial components the three vectors are, 
by Equation 1,5, 

A - A z i + A y j + A; k, B ~ iy and C - C x t + CJ t 

since this choice of axes makes By ~ B £ « C £ - o, Ex¬ 
panding the vector product (B x C) by Equation 2.20, 

(B x C) - B x Cy k j 
and also by the same rule, 

A x (B x C) - A u B x Cy i - A x B x Cyl 
Out from Equation (2.10) 

A C - A Z C Z + A M C y and A B - A z B zt 

so that 

A x (B x C) - B - A z B x Cy\ 

- (AyC tJ + A z C x B z i - A z B z C v l 

- (A Z C Z + A y Cy) B - A z B x [C z i + CJ). 

Substituting from above, 

A x (B x C) - (A‘C)B - (A B)C. . . (2.22) 
Each term of the product involves the external factor A 
in a scalar product, fust with the extreme and then with 
the middle factor. Hence the value of a vector triple 
product is entirely different if the order of the factors be 
interchanged. 

Products of more than three vectors do not often occur 
in physical applications and when encountered they are 
easily reduced by use of the preceding theorems. The 
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results serve to confirm what is already shown by Equations 
(2*21), that vectors do not follow the law of association 
for multiplication, p. 3* 

8. Summary. Chapters I and II give all the essential 
rules for the manipulation of vector algebra. By extending 
the meaning of the terms addition and subtraction to the 
geometrical processes explained on p. 6, sums and differ¬ 
ences of vectors follow the ordinary laws of commutation 
and association, as shown by Equations 1*2 and 1,3. It is, 
therefore, immaterial what order or manner of grouping 
of the vectors is followed when adding or subtracting them. 
Vectors and scalars are so far in agreement, but important 
differences appear when products are considered. 

The scalar product of two vectors is commutative and 
distributive, Equations 2,1 and 2.7, and vanishes when 
cither factor is zero, exactly as in ordinary algebra. It 
differs, however, in vanishing when the vectors are per¬ 
pendicular. The vector product of two vectors differs 
considerably from the ordinary rules. It is not commuta¬ 
tive (Equation 2.13), it vanishes when the vectors are 
parallel, and it is distributive only if the order of the 
factors Is maintained unchanged (Equation 2.18). For 
more than two vectors the non-commutative nature of the 
vector product causes the ordinary law of association to be 
not obeyed, as has been demonstrated in Section 7 of this 
chapter. 

Mathematically, vectors can be defined as quantities 
that require for their manipulation a non-commutative 
algebra which has the agreements with, and differences 
from, ordinary scalar algebra that have been summarized 
above. Vector algebra is only one of many non-commu¬ 
tative algebras known to mathematicians, but it is among 
those that have the greatest practical utility. 


CHAPTER III 

THE DIFFERENTIATION OF VECTORS 


1 . Scalar Differentiation, Let V be a vector which 
is undergoing a continuous change of magnitude and 
direction. In Fig, 15 < 5 V is a small increment which gives 



Fio. 15,*—Scalar Differentiation of a Vector 


the new value V + 6 V for the vector. Suppose that V 
is a vector function of a scalar variable t ; then when t 
changes from t to l -f- fii t V becomes V + dV. The ratio 
6 V/ 5 t is the average rate of change of V with t r and as 
fit becomes vanishingly small the ratio attains a limiting 
value which is the rate of increase of V, i,e, 

SV 
fit 


] 


as St becomes zero. 


dV T . 

A - Lun 


[ 
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This is the derivative of V with respect to the scalar 
variable /. As V varies, its extremity moves over a curve 
and the derivative is a vector in the direction of the tangent 
to this curve at each point* Second, third and higher 
derivatives are obtained by similar arguments, by analogy 
with the ordinary ideas of scalar calculus. 

When V is expressed in rectangular co-ordinates, 

V - VJ + VJ + F a k ; 

V Tt Vy* V t are now functions of t. Since i, j, k are 
constant vectors, it follows that 


dt 


dV*. dV % 


-+ 

dt 


dt 





( 3 -i) 


and similarly for higher derivatives* 

A most important physical instance is when V is a 
vector of displacement and t is time ; the extremity of V 
is then the path of a moving body* The derivative 
dW/dt is the velocity along the path at any instant, and 
d 2 Yfdt 1 is the corresponding acceleration. 

2 . Differentiation of Sums and Products. If V 
is the sum of two vectors A and B, both of which are 
functions of /, then a change of t to t 4 - dt makes 

V + 6 \ - (A + dA) + (B + dB), 
so that dV » dA + dB, 


Dividing by 6 i and proceeding to the limit gives 


dV 

dt 


-a< A + B >- 


dA d& 
dt + dt ’ 


■ • (3-2) 


so that the operation of differentiation is distributive, as 
in ordinary calculus. 

For the scalar product V - A-B the increment in t gives 
V+$V=(A+MHB+< 5 B)-A-B+ 5 A*B+A-<SB+( 5 A-< 5 B, 
expanding the scalar product by Equation 2.8. Subtract- 


THE DIFFERENTIATION OF VECTORS 31 

ing A-B, dividing by dt and proceeding to the limit, 
neglecting second-order terms, 



dA 

dt 


B + A- 


d\\ 

df 


(3-3) 


The order of the factors can be changed at will, i.e. the 
operation is commutative* If B - A we have the useful 
form for the self product, 


J)(A-A) - J t i^) - j t {A*) - 2A-~ - zA~ . (3.3a) 


dt 


For the particular case of a constant vector A* = A* — constant, 
defining the position of a point on the surface of a sphere of radius 
A, dA/dt is zero and Equation 3.3*1 leads to A*(dA/dt) «• o* 
Hence dA/dt is perpendicular to A. If a point moves on the 
surface of a sphere its velocity is always normal to the radius 
vector, as is physically obvious* 


I 1 or the vector product Y » A x B expansion by Equa¬ 
tion 2*19 gives 

V + < 5 V - {A + ( 5 A) X (B + < 5 B) - A x B + dA x B 

f Ax < 5 B + dA x dB, 

leading to 


Tj\ r> 4 dB 

Tl (A X B) = — X B + A X 


dt' 


(34) 


in which the order of the factors must be strictly maintained* 
Putting p = djdt the reader can easily verify that for 
triple products 

p[A|B x C)] - /)A-(B x C) + A(pB « C) + A(B x pC) f 
and 


^[A x (B x C)] -pA x (B x C) + A x f pB x C) 

+ A x (B x pC). 

In general, differentiation in vector analysis is seen to 
follow the same rules as in ordinary differential calculus, 
except so far as vector algebra differs from scalar algebra 
m the non-commutative property of a vector product* 
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3 . Partial Differentiation* These simple properties 
of differentiation as applied to vectors can be extended to 
partial derivatives when a vector is a function of more 
than one scalar independent variable. The most useful 
case is that of a vector V which is a function of the car¬ 
tesian co-ordinates a?, y, 2 of a point in space. If y and 2 
remain constant while x increases, the partial derivative 
dV/dx denotes the rate of increase of V with respect to x . 
Likewise changing y and 2 alone gives the partial deriva¬ 
tives dV/dy and dY/dz t denoting the rates of increase with 
respect to y and z respectively. If now x t y and z change 
simultaneously by differential increments dx, dy y dz , the 
total change or total differential of V will be 

-v av, av. 3v , , 

dv ~te dx+ d/ y +te d *’ ■ • (3 ' s) 
which is of frequent occurrence in physical applications of 
vector analysis. If r - jrf + y} + zk is the radius vector 
from the origin, then its differential increment is 


dr - dx i + dyl + 


Equation 3.5 may be symbolically written as 

** ■ [si 1 ** * + 5 ? fa ] v ■ 


If wc now define the operator V by 


V 


.3 .3 

“ 3 * + hy 



it is easy to verify that the scalar product of V, regarded 
as a kind of vector, with dr gives the operator in square 
brackets. Thus, 

rfV-(V’*)V. . . . - (3*s m 

The operator V is of immense importance in physical 
applications of vector analysis, wherein it appears in 
association with both scalar and vector operands. These 
uses will be fully explained in Chapters XV and V. 


CHAPTER IV 


THE OPERATOR V AND ITS USES 

1, The Operator V. The differential operator V was 
introduced by Sir William Rowan Hamilton and developed 
by P* G. Tait; it is of centra! importance in all three- 
dimensional physical problems. The symbol V was origin¬ 
ally named 1 nabla * after a harp-like ancient Assyrian 
musical instrument of similar shape ; other writers have 
called it 4 atled \ i.e. 1 delta ’ reversed. It is now usual 
to adopt the term * del 1 introduced by ]. Willard Gibbs. 
In cartesian notation 

^ - V - 4 + 4 + 4 ■ • • ( 4-0 

which may be applied as a directive differentiator either 
to a scalar or to a vector function of space. Again, treat¬ 
ing the differentiators in V as scalars, we may formally 
regard V as a vector which can have either a scalar or a 
vector product with other vectors. In vector analysis 
there are three fundamental operations with V which are 
of physical interest. If S is a scalar function and V a 
vector function of space, these operations are (i) V^>, where 

V acts as a differentiator ; (ii) V-Y, and (iii) V x V, where 

V is treated as a formal vector. (Sec also p. 113}. 

2 , The Gradient of a Scalar Field* On p. 12 it 
has been shown that certain physical quantities, such as 
temperature or electric potential or any such non-directed 
quantity, can be represented from point to point in space 
by a scalar point-function, S t of the co-ordinates* The 

33 
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entire scalar field can be mapped out by level surfaces* 
upon each of which the scalar function S has a constant 
value* Consider two such surfaces very close together 
and examine a small portion of them in the neighbourhood 
of a given point A on the surface characterized by the 
constant value S of the scalar function ; the second surface 
is specified by a constant value S + dS. This is shown, 


unit 
normdJ 

n 



Fig, i 6.—Gradient of a Scalar Point Function 

much magnified, by Fig. 16* If r is the radius vector 
from the origin to the point A t any point such as B in the 
second surface is given by r + dr. The least distance 
between the surfaces will be AC f in the direction of the 
unit normal vector n at A and of length dn. 

If dr be the length of AB t the magnitude of the rate of 
increase at A of S in the direction of AB will be dS/dr 
when the two surfaces are vanishingly close together* 
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This rate of increase becomes greatest in the direction of 
the unit normal n, Le. along AC t when it has the value 
dS/dn. Note that 


as 

dr 


=- cos 0 * 


Hence, if n is the unit vector normal to the level surface 
at any point in a scalar field, the vector n dS/dn gives the 
greatest rate of increase of S at the point in magnitude 
and direction* This vector is called the gradient of S at 
the point and is written 

grad S m .(4,2) 


Thus* the gradient of a scalar field is a vector field, the 
vector at any point having a magnitude equal to the most 
rapid rate of increase of S at the point and in the direction 
of this fastest rate of increase, Le* perpendicular to the 
level surface at the point* 

A simple physical example will fix the reader’s ideas* 
Suppose 5 is the potential in an electric field due to static 
charges* Then the electric force at any point is in the 
direction of the greatest rate of decrease of potential, Le* 
normal to the equipotential surfaces, snd has a magnitude 
equal to that rate of decrease. That is, the electric force 
Is - grad S and is obviously a vector field* 

It is clear from the way that the idea of gradient has 
been introduced that it is an intrinsic property of a scalar 
field and is, therefore, a physical notion entirely inde¬ 
pendent of any particular system of co-ordinate axes* 
In other words, the operation denoted by grad is invariant, 

2 a. The Operation VS. Consider now the vector 
represented by 


_ „ dS. dS t dS t 
- S 1 + Ty > + ai k 

in rectangular co-ordinates. The vector rates of increase 
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of S in the directions of the axes of *, y and 3 are 135/9x, 
idS/dy and k35/3;sr* Their sum will be a vector with 
the magnitude and direction of the total or most rapid 
rate of increase of S. It remains to show that this expres¬ 
sion is equivalent to the gradient as just defined. To do 
this, take the scalar product on both sides of the gradient 
equation with an element of radius vector dr t Le. resolve 
dr in the direction of the normal at A in Fig, 16 , Then 

(grad 5)-dr - - g^dr cos 0 - an - dS t 


since 35 / 3 n is the total normal rate of change of 5. 
in rectangular co-ordinates. 


dS 


35, 35, 35, 

dx + =-dv + ^~ds ; 
dx dy J ds 


Now 


hence, 


, _ _ . 35, 35, t 35. 

(grad S)’dr - ^dx + ^dy + -^dx 


- (‘fi + |?] + |fkY(Ad + dyi + eh k) - (VSJvfr; 


K dx 
bo that 


J „ _ _ 3 S. 5S, 3S, 

grad S - VS - . 


(4.3) 


showing that the operations grad and V applied to a 
scalar point-function are identical. Also since the gradient 
in a scalar field is independent of the nature of co-ordinates, 
so also is the operator V ; it is a mere convenience, there¬ 
fore, to express it in terms of x t y and z. We may regard 
V as a directional differentiator which, when applied to a 
scalar function of space, gives the greatest rate of change 
of the function in magnitude and direction at every point, 
i*e. derives from a scalar field its vector field of gradient* 
We have here an important case in which a vector field 
is derived from a scalar field by the process of finding 
the gradient of the latter* It does not necessarily follow 
conversely that all vector fields can be expressed as the 
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gradient of a scalar function* Let V 5 be a vector which 
is derived from a scalar 5 in the form 

Y s - grad 5 - V 5* 

In Fig* 17 draw any path, such as that marked Path i f 

V t =grad S 



FlO. 17*—Line Integrals in a Lamellar Fietd 

between two points A and B in the vector field and let 
Vj make an angle 6 with the element d\ of the path ; 
then the product of the length of the element and the 
component of V# in its direction is, from p* 18 , 

V s cos 0 dl - Vg’dl. 

The path is traced out by the extremity of a radius vector 
r from the origin ; then dl is equivalent to dr and we may 
write, from p* 36 , 

Ys'd I * (V&) dr = (grad S) dr - dS * 

Hence the line integral of V s from A to B is 
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where S A and S B are the values of the scalar at the ex¬ 
tremities of the path. Since only the end values are 
concerned, the line integral has the same value inde¬ 
pendently of the path followed. 

Now take a closed path, consisting of Path i from 
A to B followed by Path 2 from B to A. From Equation 
( 4 , 4 ), since the limits of integration are reversed for the 
second path, the value of the line integral along it is 
S d - S B ; hence round a closed path 

<j>(grad = o . . ♦ * ( 4 . 5 ) 

Summarizing, when a vector field can be expressed as 
the gradient of a scalar field, the line integral of the vector 
taken between two points is independent of the path 
followed and is equal to the difference between the values 
of the scalar at its ends ; further, the line integral round 
any closed path in such a vector field is zero, 

A vector field V 5 derived from a scalar S by the relation 
Y s = grad S is sometimes called a scalar potential field , 
5 being the potential of Since space is divided up 
into layers or laminae by the level or equipotential surfaces 
of the function S , V 5 is also called a lamellar vector , 
Again, since the essential property of such a vector is 
that its line integral round any dosed path is zero, the 
vector field is called a non-curl field , a term which will 
be explained in Section 4 . 

A simple example is the electrostatic field of charged 
conductors, which can be specified by a system of equi¬ 
potential surfaces. The vector field of electric force is 
derived from the scalar potential as its negative gradient, 
i.e. the electric force is in the direction of the greatest rate 
of fall of potential and has a magnitude equal to that rate. 
The line integral of the electric force between two points 
is the potential difference between them and is independent 
of the path taken ; it represents the work done in moving 
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a unit charge of positive electricity from one point to the 
other. Round a dosed path the work done is zero, ex¬ 
pressed by the statement that the electromotive force 
vanishes round any closed path in a static electric field. 
The reader will be able to construct other physical examples 
of lamellar vector fields, e.g. in the flow of heat with tem¬ 
perature isothermals or in gravitational attraction with 
level surfaces. 

3. The Divergence of a Vector Field. In Fig. 18 



Fig. 18,—Divergence of a Vector Point Function 


let V be the value of a vector function at the middle of 
an infinitesimal element of volume with sides dx t dy and 
dz parallel to the axes of jc, y and z* The vector V has 
axial components of magnitude V y and V z . To fix 
ideas the reader may think of V as a vector giving the 
velocity of a moving fluid in magnitude and direction. 

Consider the two faces of the volume element, each 
with area dy dz perpendicular to the axis of X, On the 
left-hand face the value of the x component of the vector 
at the middle of the face becomes 


Vi ~ hx dx> 


4 
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and may be taken as the value all over the face when this 
becomes vanishingly small. Similarly, on the right-hand 
face the x component is 


11 + *a****■ 


We now define the flux through any face as the scalar 
product of the vector area of the face and the vector V ; 
he, as the product of the area of the face and the normal 
component of the vector upon it, see p. ig. The flux 
is positive when the component of V and the outward- 
drawn normal on the face of the element are in the same 
sense. Then the excess of flux leaving the dement over 
that entering it in the X direction is 

^Vz+^-j^dx^dy dz-^V^^~dxjdy dse***~dx dy dz. 

In the hydromechanical case this represents the net volume 
of fluid passing per second in the X direction. By similar 
reasoning the contributions parallel to Y and Z are 

-jfyd* dy dz and ^-^dx dy dz . 

The total net flux diverging from or leaving the dement 
is, by Equation 2 . 12 , 


/m 

\dx 


+ 


dy dz 


^jdx dy dz. 


The amount of this flux per unit volume is defined as the 
divergence of the vector V and is written 


div V- 


3T* dV v 

dx + dy 



(4-6) 


Since the divergence is the amount of flux, it is essentially 
scalar. 

If the divergence exists at a point in a fluid, whether 
liquid or gas, and is positive, it expresses the rate at which 
fluid is flowing away from the point per unit volume thereat. 
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Hence, either the fluid is expanding and its density at the 
point is falling with time, or the point is a source at which 
fluid is entering the field. When the divergence is nega- 
live it gives the rate at which fluid is flowing toward the 
point per unit volume. In this case either the fluid is 
contracting and its density rising at the point, or the 
point is a negative source, a sink , at which fluid is leaving 
the field. Since most practical liquids are almost incom¬ 
pressible, the existence of divergence in them means the 
presence of a volume distribution of sources or sinks rather 
than changes of density. In the case of non-material 
fluxes, such as those of the thermal, electric or magnetic 
fields, the existence of divergence means the presence of 
a source or sink of flux at the point. For example, in the 
electric field, positive divergence means that there is 
positive electricity at the point; in the thermal field the 
point is either a source of heat or a place where the tem¬ 
perature is falling. 

When the divergence is everywhere zero, the flux enter¬ 
ing any element of space is exactly balanced by that leaving 
it and %ve may write, 

div V - 0 ,.( 4 . 7 ) 

which is true in many practical problems. In a fluid 
this means that there can be no sources or sinks in the 
field, nor can its density be changing ; i.e, the fluid is 
incompressible. If the fluxes entering and leaving an 
clement are equal, none can have been generated within 
it; the lines of flow of the vector V must either form 
closed curves (cf. the magnetic field of a current), or 
terminate upon bounding surfaces (cf, the electric field 
m a condenser), or extend to infinity. A vector which 
satisfies this condition is allied sohnoidal (from a Greek 
word meaning a tube). 

3£i* The Operation V*V. Consider now the scalar 
product of the operator V aitd the vector V. Expressed 
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in rectangular co-ordinates the rules for scalar products, 
Equations 2,5, 2.6 and 2.$ give 


v-v 


■( i i +, ! + k s )' ( ’ y+i ' j+r,k)s5 


d I*JlyJX± 

dx + dy ox ’ 


that is, 

V ‘V = div V .... {4.8) 
Now the operation V has been shown on physical grounds 
to be invariant, i.e. independent of any system of axes* 
Hence the idea of divergence is also invariant* This is, 
indeed, physically obvious, since the amount of flux enter¬ 
ing or leaving per unit volume at a point is clearly a con¬ 
ception quite independent of any system of co-ordinates 
that may be used to express the shape and position of 
the volume dement. 

4 * The Curl of a Vector Field. It has been shown 
on p. 38 that when a vector field can be derived as the 
gradient of a scalar field, the line integral of the vector 
taken round any closed path is zero ; this result is true 
no matter what size or shape the path may have. A 
vector field satisfying this condition is known as a lamellar 
held, and is of a special, though very important, class. 
Many vector fields occur in physical problems, however, 
in which the closed-path line integral is not zero and 
which cannot, therefore, be expressed as the gradient of 
a scalar point function ; it is to an important property 
of these more general fields that we now give attention. 

Consider a very small region of such a vector field, 
several lines of flow in which are shown by Fig. 19(d) ; 
the portion is chosen small enough for the lines to be 
regarded as nearly straight and parallel. Into this field 
put a small plane area, shown for convenience as a rect¬ 
angle. When the area is perpendicular to the field, as 
in position 1, none of the field is directed along any part 
of the bounding edge of the area ; the line integral round 
it is zero. In position z, with the area parallel to the 
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field, since the value of the vector along the upper edge 
is assumed to be different from that along the lower edge, 
the line integral round the boundary has a finite value. 
Similar arguments apply to intermediate positions ; the 
value of the line integral depends, therefore, upon the 
direction of the normal to the area relative to the field. 



Le. upon the orientation of the given small vector area at 
the region considered. 

In general, if we put a small vector area of any shape 
at any point in a vector field and compute the line integral 
of the vector V around its bounding edge there will be 
an orientation of the area for which the line integral is 
greatest. The amount of this maximum line integral ex¬ 
pressed per unit area is called the curl of the vector field 
at the point, and is given the vectorial sense of the positive 
normal drawn on the small exploring area when in the 
position giving this greatest Integral. Some writers use 
the term rotation (briefly rot.) to denote the same con¬ 
ception, since the curl is associated in hydromechanics 
w ith the rotation or spin of a fluid in which the particles 
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have angular velocity. Physical examples illustrating the 
meaning of the curl of a vector held will be given later ; 
the meaning depends very much upon the nature of the 
problem. 

To calculate the curl in terms of its cartesian com¬ 
ponents, take three infinitesimal rectangular areas inter¬ 
secting mutually at right angles at a point where the 
vector field V has components of magnitude V It V v , v„ 
as in Fig. ig(6), Taking the positive normals to the areas 
along the positive directions of the X f Y and Z axes 
respectively, the circular arrows indicate the positive senses 
in which their boundaries must be traversed to accord 
with the right-hand screw rule for vector areas. 

As an example, take one of these areas such as abed in 
Fig. I9(fc) with sides dx , dy\ its normal being along the 
axis of Z, Since the rectangle is very small, the numerical 
value of the component of V at the middle of any side may 
reasonably be taken as the average value along that side; 
the arrows show the directions in which the components 
act. Since V xt V !/f V x are functions of the co-ordinates 
(#, y i ~) of the middle of the rectangle, the average values 
along the four sides ab t be , de t ad will be respectively 
-3 Vx, Tr rSVy dVtj , T7 t dV u , 

v * - i-dy d y> v » + i ^ v ** i-fa d y and V v- tr dx dx - 

Around the contour abed the line integral is, therefore, 

that is* 

Since the area of the element is dx dy , the bracketed term 
is the magnitude of the component curl of the vector field 
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taken about the Z axis. Giving it the sense of unit vector 
k, therefore, we may write 


curl, V 



By following a precisely similar method with the two 
remaining rectangles we find the component curls about 
the Y and X axes to be 


curl r V - and curlj V 



Adding the three components gives 

curlV-ffi dV Au( dV * dV *\lJ dV u dV *V / x 


which can be conveniently written in the form of a deter¬ 
minant as 


curl V 


I I k 

d/3x d/dy d/dz . 

K v v v. 


( 4-9 s “) 


4 a. The Operation V x V. Take the vector product 
\ r x V of the operator V and the vector V, expressing 
the result in rectangular co-ordinates. Using Equation 
2.2o for expanding a vector product, 


V * V ‘ ('si + >fy * k I) « (W + + »W 



that is. 


V x V - curl V . , . (4*10) 

"I he operation is, of course, independent of the system 

of axes, 

5 , Simple Examples of Curl. The meaning of the 
operation curl may now be illustrated by a few simple 
examples. First consider a rigid body which is rotating 
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with an angular velocity co about an axis OA , O being a 
fixed point in the body. Then any point P, Fig, 20, 



moves in a circular path about OA with a tangential linear 
velocity cor sin 0, where r is the distance from O to P . 
In vector notation the angular velocity is a vector to 
drawn along the axis in a sense related to the rotation by 
the right-hand screw rule ; the co-ordinate of P is the 
radius vector i\ Then the tangential velocity at P is 
oj x r, which has magnitude cor sin 0 and is perpendicular 
to the plane containing to and r* If in addition the 

* To verify that it is physically correct to treat an angular velocity 
as a vector it is only needful to show that two such velocities 
compound as vectors, see p. 7* Let to* be the angular velocity 
about an axis through G ; then the tangential velocity at P h 
to, x r. Now let to, be the angular velocity about a second 
axis through O inclined to the first; then Phas a tangential velocity 
due to this rotation of to, x r. But linear velocities compound 
vectorially, hence the total linear velocity at P is 

Wi xr + w,xr - (to t 4 o>,) x r = <0 x r, 
by the distributive law, Equation 2,1$* Hence the motion is the 
same as that due to an angular velocity equal to the vector sum of 
the two components. A similar theorem is obviously true of 
infinitesimal angular rotations and of angular accelerations ; it 
is not true of finite angular rotations. 
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whole body is moving with a linear velocity v 0 in any 
direction* the total velocity at P is 

V - v 0 + to x r, 

and curl V - V x V ~ V x v 0 + y * to x r. 

Now since v 0 is constant for all points in the body and is 
thus independent of x t y s z, V x v 0 is clearly zero. 

I he angular velocity to is also a constant vector for all 
points and can be written as 

to - q) x I + coj + m 2 k, 

where its components co xi co z arc independent, there¬ 
fore, of the coordinate r of the point P, Writing 

r - xi + yi + zk t 
and using Equation 2.20 makes 
to x r - {m y z - co 2 y)i + (<o z x - oi x z)i 4 {co^y * cj y x)k. 

Again, using Equation 2.20 with Equation 4.1 and remem¬ 
bering that since co xt co y% co z are not functions x,y y z their 
derivatives are zero, it is easy to verify that 

V x (to X r) - 20) x i + 20)J + - zu>. 

Finally, therefore, 

curl V - 2to 

and OJ - 4 curl V.{4.11) 

Hence when a rigid body is in motion the curl of its 
linear velocity at any poiot gives twice its angular velocity 
in magnitude and direction. (See also p. 113.) 

Xow examine the motion of a fluid and consider what 
may happen to an infinitesimal volume within it. The 
volume can have three kinds of motion simultaneously : 
(0 It may be moving with a linear velocity of translation 
23 a whole, (ii) If the fluid is deformable the element 
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may change its shape, (iff) It may be in rotation. At 
any instant the little volume may be regarded as a rigid 
body; the curl of the velocity of the fluid at the point 
where the element is situated gives twice its angular 
velocity, the instantaneous axis of rotation being that of 
the curl. The nature of the rotational motion will be 
made clearer by Fig. 21(0), which shows two positions of 


(«) 0) 

Fic. 21.—Rotational (a) and Irrotational (6) Motion 

a small portion of the fluid in movement about an axis 
at O, It is clear that the portion has rotated ; and if 
every elementary portion of the volume immediately round 
O has rotated by the same amount, curl V would give 
twice the angular velocity of rotation about 0 . By con- 
trast. Fig, 21(6) shows a small portion of fluid which, in 
its motion about O, docs not rotate ; hence there is no 
curl V and its angular velocity is zero. 

If a motion is such that the velocity has a curl, the 
ultimate particles of the body are in rotation with an 
instantaneous angular velocity. The motion is described 
as rotational or vorticaL If, on the other hand, the par¬ 
ticles do not rotate, there is no curl and the motion is 
irrotational or non-vortical. 

Another example of an entirely different kind occurs in 
the magnetic field of a conductor carrying a steady current* 
At any point in the field put a very small plane area and 
turn it into such a position that the line integral of the 
magnetic force taken round its boundary is the greatest 
possible ; this value expressed per unit area is the vector 
curl H, he, the magneto-motive force per unit area at 
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the point. If the point is within the material of the 
conductor at a place where the vector of current density 
is c, then this will be the total current passing normally 
per unit area through the elementary path when the line 
integral round it is greatest. In electromagnetic units, 
curl H — qjte r so that the curl of the magnetic force at 
any point is proportional to the current density thereat 
and has the same direction. For a point in the field 
external to the conductor there is no current density and 
curl II - o. 


6 Divergence of a Vector Product* I n certain physical 
problems, e.g. in calculating the flux of energy in an electromag¬ 
netic field, it is required to find the vector product of two vectors 
and then to work out the divergence of the result. Treating this 
at present as an exercise in manipulation of vector operations, we 
are to find div (A x B)* Using Equation z.zo, write 

V - AxB- AiB#) i+( A-B x -A x B s )j + (A x B v -A M B z )k 
-F a i+F l J + F I k* * 

Now apply Equation 4,6 and 

di vv,*^ + ?£*+£*:* 

dx +dy + d* 


. 4 n ( 8A * dA A 4. 7 ? BAA 

\ By dx) + M W “ di) + Hsr" -dy) 

- aJ^ z - - d -fy\ _ A ( dB * BBt \ A dB A 

By dj~ A *[ ~Bi - srj - - -ft)* 


Using Equations zao and 4,9 gives 


div {A x B) - B* curl A - A* curl B , . (^ lz y 

* * Divergence and Curl of 5 A. In later work we require 
the divergence and curl of a vector field which is itself the product 
of a scalar field S and a vector field A at every point. The com¬ 
ponents of the product are SA X> SA U and SA ft 
Using Equation 4.6 to find the divergence, 

div SA - JW) + f y (SA v ) + l(SA t ) 


_ s( d -i-+ + dA ‘\ . ( A bs ds ds\ 

\ dx + dy + -3z) + [ A 'dx + + A ‘Tz) 

™ S div A + (AJ + AJ + Ak){ffi + gj + gk). 
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which from Equation 4.3 is 

div SA - S div A + A* grad S , 1 (413) 

Using Equation 4.9 bis the curl can be written 

l i I * 

curl SA — | d/dx 8 /dy B/Bx , 

j SAx SAy SA Z 

the x component of which has the magnitude 

Hy<™ - | <M >> ‘ S (w-’-£) + [ A f " - 1 '!]' 

Using Equation 4.9, the first bracket is the x component of curl A \ 
the second bracket is the x component of (grad S} x A, as may be 
verified from Equations 4,3 and 2.20. Working out the other 
components gives finally 

curl SA — S curl A + (grad S) x A 

& S curl A-Ax grad S . * (4.14) 


CHAPTER V 

FURTHER APPLICATIONS OF THE OPERATOR V 


1* The Operator div grad- In Equations 4,2 and 
4.3 it has been shown that if S is a scalar function of 
position in space, then grad S==yS is a vector function 
having the magnitude and direction of the greatest rate of 
increase of S; this vector is the gradient of S and the 
flux lines of grad S cut normally through the level surfaces 
of S t Since grad S is a vector it can have a divergence ; 
using Equations 4.8 and 24 this scalar quantity is 

div grad S - V'(VS) - VV(S) - V 2 £> 

From Equations 4.1 and 4.3, 


r a e A 3 1 9 s d \f SS * BS t dS t \ 

d!V grad S = (i Fx + + k g -ji + g-i + 9jkj. 

Expanding by the usual rules for scalar products, Equations 
2.5, 2.6 and 2*8, 

div grad + + . . . (s,) 

and the operation 

3 s d 2 d % 

divgrad =^ + -5 + ^-V 2 , • . (5-2) 


which is known as Laplace's operator. Examples of its 
use will be found in later chapters. 

2 * The Operator curl grad. Since grad S is a vector, 
it is possible to calculate its curl and thus interpret the 
operation curl grad applied to a scalar point-function. 
Now from Equation 4,9 bis and above* since the com- 

5 1 
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ponents of the vector for which the curl is to be found 
are dS/dx, dS/dy and dS/dz, 

i j k 

curl grad S d/dx d/dy d/dz = zero * {5-3) 

-V *(VS)~ dS/dx 3 S/by dS/dz 

This result follows at once from physical considerations. 
On p. 38 it has been shown that = grad S is a lamellar 
vector, of which S is the scalar potential. The charac¬ 
teristic feature of a lamellar field is that the line integral 
of Vs round ajiy closed path is zero. Since the curl of 
a vector field is a particular kind of closed-path line 
integral, curl V^ - curl grad S is necessarily zero as V$ 
is lamellar. A typical non-curl field of this class is the 
electric force due to static electric charges. 

3 . The Operator V a with Vector Operand. A 
vector point-function V may be expressed in terms of its 
rectangular components F^i, V y }, F £ k in the usual way. 
Since V xr V y and F* are each scalar functions of position 
(x t y\ 2)1 Equation 5.2 applies to them all. Then if 

V-Fl + Ff + Fk 

v*v - vV*t + v%]+ v*vx . . (54) 

which is a result of great importance in electromagnetic 
theory and hydromechanics, 

4 . The Operator grad div* If V is a vector field, 
div V is a scalar field, which, therefore, has a gradient. 
This new vector, of which div V is the potential, is neces¬ 
sarily lamellar because curl grad on scalar operand is 
zero ; it may be written as 

grad div V~V(V * V )-(4 + i| +k ^)(^ + ^ + ®) 

( dw x d a ~V v d 2 V t \. (dW x dW y 5 3 FA. 
y 3 a’ 2 + dxdy + dxdz J 1+ ^ 9 # 9 y dy* dydzj 

<*> 
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in cartesians. For example, if V is a fluid velocity, div V 
gives the scalar rate at which the density at a point is 
changing per second ; grad div V gives the magnitude 
and direction in space of the greatest rate of space increase 
of the density. In an electric field V is the electric force, 
div V the density of space charge and grad div V the 
greatest rate of increase of charge at any point. 

5 . The Operator div curl. If V is any vector field, 
curl \ is also a vector field, for which we may calculate 
the divergence and the curl. Dealing with the first of 
these, using Equations 4.6 and 4.9, 



Hence the curl of any vector field is itself a solenoidal 
field, i.e, one in which the lines of flow form closed curves, 
since it has no divergence (see p. 41). 


6. The Operator curl curl. The y and * components 
of curl V are 


dz 


_dVA 

dx j 


J and 


'dV_ v 

dx 


dV, 

(£y 


>* 


and hence by Equation 4.9 the a; component of curl curl V 
will be 


curb curl V 



Expanding the factors, adding and subtracting dWJdx* 
gives 


curl* curl V 


JW£\ i_r iv 5*., 3H r/ . 

dx* dxdy eirSeJ |_3* !S+ 3y s+ 3!**J 


which from Equations 5.5 and 5.4 may be identified as 
die * components of grad div V and V 2 V respectively. 
Permuting i, j, k; x, y, z; and V x , V y , V, in cyclic 
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order gives curl, curl V and curl, curl V respectively; 
adding the three components 

curl curl V - V x {V x V) = grad div V - V 2 V, (S- 7 ) 

so that the operation 

curl curl = grad div — VN * ■ • ( 5 '®) 

a theorem which is much used in various problems. 

7 . The Classification of Vector Fields. We have 
distinguished so far between only two kinds of vector 
fields. The first is characterized by zero curl and can 
be expressed as the gradient of a scalar potential, p. 38 ; 
it is, therefore, known as a lamellar or scalar potential 
field! The second has zero divergence; there are no 
sources of vector flux in the field, the flux lines arc closed 
curves and the field is termed solenoidal, p. 41. Since 
the operation div curl on a vector operand yields zero, 
Equation 3.6, a solenoidal vector can always be expressed 
as the curl of another vector ; see below. _ 

In practice there are four kinds of vector fields, having 
different associations of curl and divergence ; they may 

be classified as follows: _ , 

(i) curl V - o and div V - o : Since the field is lrrota- 
tional or lamellar the first condition makes V - grad S, 
The second condition makes the field also solenoidal ; 
from Equation 5.1, div V = div grad 5 - V 2 S= °> which 
is Laplace's equation. Typical fields of this sort are the 
electric force in free space due to static charges on boundary 
conductors, or the irrotational motion of an incompressible 
fluid. See Chapter VIL 

(ii) curl V - o but div V * o : Here the field is still 
irrotational but is not solenoidal. Hence V *S * o, which 
is Poisson's equation. The field must contain a spatial 
distribution of sources of flux. Examples are the electric 
field of a volume distribution of charges, e.g. electrons in 
a thermionic tube, or the gravitational force inside a mass. 
See Chapter VII. 
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(m) curl V + o but div V = o. The field is rotational 
and cannot have a scalar potential, but it is solenoidal 
Let A be a vector such that V « curl A; then from 
Equation 5.6, div V = div curl A = o as required. Equation 
5,7 makes curl V = curl curl A = grad div A - V B A. The 
nature of the vector A is a matter of choice \ suppose, 
therefore, that we desire it to be solenoidal, then div A <* o 
and 

curl V * - V 2 A + o* 

This is similar in form to Poisson’s equation but contains 
a vector instead of a scalar operand. For this reason A 
is termed the vector potential of V. Important typical 
cases are the magnetic field within a conductor carrying 
a steady current and the rotational motion of an incom¬ 
pressible fluid. See Chapter VIII. 

(iv) curl V * o and div V + o. This is the most 
general type of vector field possible, found in the rotational 
motion of compressible fluid and in some branches of 
electromagnetism. The field can be expressed as the sum 
of two fields ; the first is a lamellar vector, which has no 
curl but can have divergence, and the second is a solenoidal 
vector, which has no divergence but can have curl. Writing 
V - grad S + curl A, 
div V = div grad S = 4- o, 

from Equations 5.1 and 5,6. This is Poisson’s equation 
and determines S. Using Equations 5.3 and 5.8 and 
choosing A to be solenoidal, as before, 

curl V = cur! curl A = - V 2 A + o, 

which determines A, Thus a general vector field may be 
decomposed into two fields, one being lamellar, with 
scalar potential S t and the other solenoidal, with vector 
potential A. This is sometimes known as Helmholtz's 
theorem , 

i 32 illustrates diagrammatieally the characteristics of vector 
odds near a point. Supposing the vectors to illustrate th^ 
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morion of fluid particles, Fig. 22(a) shows neither divergence nor 
curl. Fig. 22(b) shows divergence alone ; fluid is streaming from 
the point, at which the density is falling ; alternatively the point 
is a source. Fig. 22(f) shows fluid in rotation about the point ; 
there is curl but no divergence, and the fluid is incompressible. 
Finally, Fig. 22(d) illustrates particles in a compressible fluid in 
rotation ; there is both curl and divergence. 



Q />/ 


no divergence divergence curt 
no curl only only 


divergence 
end cur / 


M (4) M (<9 

Fig. 22.—Typical Characteristics of Vector Fields 


8, Two Useful Divergence Formulae* In the application* 
of vector analysis two special forms of the divergence formula, 
Equation 4.13, are often used. 

Let U and W be two scalar point-functions, grad W being a 
lamellar vector field with W as its potential. Then in Equation 
4.13 put U for S and grad W for A, giving 

div (U grad W) - UV'W + V £Ly W 

- Uy *W + (grad £ 7 )'(grad W) . (5.9) 

Again, let 5 be a scalar field and A = curl B a vector field, 
B being the vector potential of A ; then Equation 4.13 becomes 

div (S curl B) — S div curl B + (curl B)*(grad 5 ). 

But Equation 5,6 makes the first term zero ; hence 

div (S curl B) — (curl B) (grad S). . , ( 5 -to) 


9 , The Vector Field grad (k/r), One of the most 
important lamellar fields is that which has a scalar potential 
S = k/r, k being a numeric and r the length of the radius 
vector r - xi + y] -f tfk from the origin to any point P 
in the field. The level surfaces are concentric spheres 
with a common centre at the origin ; the potential increases 
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toward infinity as the origin is approached. From Equa¬ 
tion 4.3, 

grad S - grad (k/r) - V(*/r) 



- - k(xl 4- yi + zk)/(x 2 + y* + * 2 ) s'* 

Simplifying, 

grad (k/r) - - kr/r* - - ftr,/r\ . . ( 5 .u) 

where r, — r/r is unit vector in the direction of r. As 


spherical 



Fig. 23.—The Point Source 

shown by the minus sign the gradient is directed radially 
inwards, as drawn in Fig. 23 ; the lines of flow of the vector 
field grad (k/r) are thus radial lines cutting perpendicularly 
through the spherical level surfaces and converging upon 
the origin. Excluding for the moment the origin itself, 
at which k/r and grad (k/r) both become infinite, calculate 
the divergence at any other point by means of Equation 5.1. 
Then the reader can easily verify by direct partial differen- 
Hation of k/r that 

VW') " div grad (k/r) - o, . (5.13) 

*0 that the potential satisfies Laplace’s equation, p. 54, 
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everywhere except the origin, where the divergence is 
infinitely great. 

Hence when the scalar potential is i/r, its gradient is 
-i-i/r 2 . This is the familiar case of the unit point source 
or sink common in all branches of theoretical physics. A 
potential of this kind is known as a Newtonian potential, 
since it is that of a point unit gravitating mass at the origin ; 
the gradient gives the attractive force per unit mass placed 
at P, this force varying with the inverse square of the 
distance. Many other phenomena also follow a similar 
law. A negative electric charge at the origin attracts a 
positive charge at P; if both are of unit magnitude. 
Coulomb’s law shows Equation 5.11 to be true. Again, 
in hydromechanics the velocity of the fluid toward a point 
sink—where fluid is abstracted from the field—varies in¬ 
versely as the square of the distance and satisfies the same 
equation. 

On p. 56, if unit mass attracts a similar mass at unit 
distance with unit force, k = 1. Following the convention 
of p. 13, if a single flux line per unit area of the level 
surface is drawn perpendicular thereto to represent unit 
value of the gradient at unit distance, it follows that the 
total flux of grad (i/r) is - 4-% since unit sphere has a 
surface of 47 1 and the flux is inward. If the quantity of 
attracting matter at O is k, the flux is - 431 A, The subject 
will be further discussed in Chapter VII. 




CHAPTER VI 

THE THEOREMS OF GAUSS AND STOKES 

I. The Divergence Theorem of Gauss. Consider 
a closed surface s drawn in a vector field V ; the normal 
flux through an element of vector area ds ** nds t n being 
the outward drawn normal, is Vv/s ° V-n ds 9 as explained 
on p. 19. The surface integral of such elements will 
give the total normal flux through the surface in the form 
given by Equation 2.1 z, 

Referring to Fig. 24, an elementary volume dv within 


V 


Fig* 24.^Gauss 1 a Theorem 

the surface i Is shown ; a small cube has been taken for 
convenience, but any shape will suffice. The size of the 
element is much exaggerated, merely for ease in drawing 
the diagram. The total flux diverging from this volume 
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is div Wff, V being the vector at the middle of the element, 
and this flux has been computed on p. 40 by calculating 
the surface integral of the normal components of V 
through the surfaces bounding the volume. For the face 
abed the positive direction of the x component of V and 
the outward-drawn normal are in the same sense and the 
flux is positive, A contiguous cube, shown dotted, has 
the same component of V acting through the common 
face—since V and its components are assumed to be 
continuous, as also are the derivatives of these quantities— 
but its outward normal, shown dotted, is in the opposite 
sense and the flux is negative. The surface fluxes through 
the common face cancel. This argument can be applied 
by adding further elements along the X axis until we reach 
an element which has one face in the surface f, shown 
shaded ; this alone makes a contribution to the normal 
flux through $, Applying the same treatment to the 
surfaces of all elements of volume throughout the space 

enclosed by 5, we arrive at a total flux J j YJs, But 

at the same time we have integrated div Ydv throughout 
the enclosed volume and this also measures the total flux. 
Equating them. 


| j \-ds - J J \nds - 11 div ydv - J| | v -ydv, ( 6 . 1 ) 


which is Gauss*s theorem of divergence. 

Although the deduction of this theorem is not a rigid mathe¬ 
matical proof—many such can be found in the larger text-books— 
it k based on self-evident physical facts. Consider, for example, 
a dosed surface drawn within a fluid which is moving at a given 
point with vector velocity V, The total amount of fluid passing 
per second through the surface can be found in two equivalent 
ways. First, by calculating V*ds, i.e. the product of an dement 
of surface and the component of velocity perpendicular to it, for 
every element of the surface and adding all the contributions. 
Second, by investigating the divergence of a volume element, i.e. 
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the excess of sources of fluid over the sinks oontained in it per 
unit volume, end integrating divVdv through the enclosed 
volume, The two results are physically equivalent, since the 
excess of fluid leaving the surface over that entering it must be due 
to. the net amount of fluid injected into the portion of the field 
within i by the aggregate of sources and sinks. 

2. Gauss’s Theorem and the Inverse Square Law. 
In Fig. 85(a) let ds be an elementary portion of a surface, 
n being unit positive normal upon it. From a point O 
draw a conical pencil touching the boundary of ds and 
let r 1 be unit vector in the direction of the radius vector 
r “ n 1 from O to the clement. Draw spheres of radii r 
and unity respectively. If da is the area cut by the 
cone from the surface of the sphere of unit radius, 
rfw/i s - ds cos 6 /r* ; da is called the solid angle sub¬ 
tended by ds at O. The solid angle is regarded as positive 
when the angle 6 is acute and negative when obtuse. 
The total solid angle subtended by a spherical surface 

is 4.T. 

Fig. 25(6) shows any closed non-sphcricat surface. 
When O is inside the surface an element at 1, 2, or 3 
contributes a positive solid angle. A cone joining O to 
elements 4, 5, and 6 shows that these subtend equal solid 
angles at O, but while the contributions of 4 and 6 are 
positive that of 5 is negative. In general, therefore, any 
small cone from O cuts the surface an odd number of 
times, equivalent to a single contribution of dm. Hence 
the entire surface s subtends at 0 a total solid angle of 
When the vertex of the cone is at O' outside the 
surface, the surface is cut an even number of times with 
alternately negative and positive equal contributions, as 
shown at 7, 8, 9, 10. Hence the surface s subtends a 
solid angle of zero at an external point. These two facts 
are a purely geometric consequence of the definition of 
solid angle. 

Consider now the Newtonian potential field of a point 
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/ig, 25, —Gauss’s Theorem and the Inverse Square Law 


THEOREMS OF GAUSS AND STOKES 63 

unit mass or other similar physical entity placed either 
at O or at O', Taking its potential as 1 fr t Equation 5.11 
shows that vve can put in Equation 6.1 

V - grad (i/r) = - r^r 2 

and Y ds ^ - (r t fr 2 )-nds =■ — cos 0 dsfr z = - da) y 

noting that r* and n axe unit vectors with an angle 0 
between them. Then 

J J V.rfs - - j j dm * . * ( 6 . 2 ) 

Hence the total closed surface integral of the vector 
grad (i/r), e,g, of the total flux of force from a point unit 
mass, is - jpT when the mass is inside the surface and 
zero when outside it. The minus sign means that the 
flux is Inward, see also Fig, 23, In both cases the flux 
is independent of the position of the mass, which may be 
anywhere inside or outside the surface. This form of 
Gauss's theorem is extremely useful and will be used in 
Chapter VIL 

3 , Stokes’s Theorem, In a vector field V draw any 
unclosed surface or cap having for its bounding edge a 
given closed curve, and calculate the line integral of V 
round the closed curve when it is traced out counter- 
clockwise as shown by the boundary arrows in Fig, 26, 
Then from Equation 2,11 the value of this integral will be 

the circle indicating that a closed path has been traversed. 
By means of a system of intersecting lines drawn in 
the unclosed surface let it be divided up into infinitesimal 
surface elements. Consider the shaded element of vector 
area nds t n being the positive unit normal upon it. The 
boundary of the element is traced out counter-clockwise, 
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this sense and the direction of n following the usual right- 
hand screw rule, see p. 2 z. At the middle of ds the 
vector field will have a curl and the scalar product of this 
with the vector area will give the component of the line 
integral of V round the boundary of ds t Le. nxurl Wi, 
as explained on p, 43. A similar process can be applied 
to all surface elements, tracing them all out in the same 
way* Since V and its derivatives, and hence curl V, are 



Fig. 26. —Stokes’s Theorem 


c/osed curve 


V 


all supposed to be finite and continuous functions, the 
line integral along the common sides of contiguous elements 
will cancel, leaving only a contribution from sides that lie 
in the given closed curve. Adding the contributions of all 
elements yields the dosed line integral, but in obtaining 
it a surface integral of n curl Vds has been simultaneously 
made ; hence, 

(j> V-dl - ||n*curl Yds, # , (6.3} 

This is Stokes's theorem and states that the tangential line 
integral of a vector function round any closed curve is 
equal to the normal surface integral of the curl of that 
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function over any unclosed surface which has the curve 
for its bounding edge ; rigid mathematical proofs taking 
account of all the necessary continuity conditions would 
be out of place here and will be found in the usual text¬ 
books* The converse theorem is necessarily true, i,e. if 
the line integral of V round a closed curve is equal to the 
surface integral of A over any open surface bounded by 
the curve, then A - curl V. 

When curl Y is zero. Equation 6.3 shows that the closed- 
path line integral of V vanishes ; this has been shown on 
p* 38 to be a characteristic feature of a lamellar field* 

4 . Invariance of Divergence and Curl* Gauss's and 
Stokes's theorems may be regarded from several points of view. 
Their physical meaning has already been stressed, the first theorem 
providing alternative ways of expressing the dux of a vector through 
a closed surface, the second giving equivalent ideas for the line 
integral of a vector round a closed curve. Both results are 
necessarily scalar quantities* From another standpoint we may 
look upon them as useful analytical transformations which enable 
certain integrals to be more simply expressed* Equation 6.1 
shows that the volume integral of a vector function can be ex¬ 
pressed as a surface integral over a closed surface confining the 
volume. Equation 6.3 reduces the integral of a vector over an 
open surface to one taken round its boundary edge. 

Returning to Equation 6,t, confine the volume integral to the 
space within an infinitesimal dement dv and the surface integral 
to its enclosing surface, then 

V-ds dlv V dv , 

The divergence at a point is the value of this surface integral per 
unit volume as the element is made vanishingly small, i,e, 

v*] as dv becomes zero. 

Similarly, apply Stokes’s theorem to an infinitesimal plane area 
at a given point, turning the area until the line integral of Y round 
its boundary has the largest value ; then n and curl V are in the 


diy V — Lira, [jgj J j" 
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sunme direction and their scalar product i$ the magnitude of the 
curl, say curl V, giving 

(J)V </I = curl V du 

J mar. 

Hence, as ds is made vanishingly small, 

>fc?] 

and has the direction of n in the position of maximum line integral. 
The ideas of divergence and curl have been shown on physical 
grounds to be conceptions independent of any system of co* 
ordinates in which they may be expressed ; they are invariant. 
The formulae derived above can be regarded as a purely mathe¬ 
matical justification of this property, since no co-ordinates are 
involved in their deduction. Some mathematical treatises, indeed, 
ti eat these formulae as definitions of div and curl, building up the 
theory of vector analysis thereon. 


curl V - Lim. 


■[ 


CHAPTER VII 

THE SCALAR POTENTIAL FIELD 

1 . General Properties, Let S be a scalar point-func¬ 
tion which may be mapped out in space by a series of 
level surfaces, upon each of which the scalar has a definite 
but different constant value. These surfaces divide up 
the region of space into a series of layers or laminae. 
Associated therewith is a vector field directed every¬ 
where normal to the level surfaces, he. in the direction 
of the greatest rate of increase of S at any point and having 
a magnitude equal to that rate of increase. Tins is ex¬ 
pressed by Equation 4.3, 

Vj-gradS-VA 

Fig, 27 illustrates these conditions for small portions of 

grad S 



two infinitesimally-close level surfaces, characterized by 
values of S and S + dS respectively. 
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Referring again to Fig. 17 and Equation 4.4^ it has been 
shown that such a vector held has the essential property 
that the value of the tangential line integral of along 
any path joining two points A and B is the same and equal 
to the difference between the values of S at B and A t he* 

-S B - S A . 

Hence round any closed path the line integral of Y s is 
zero ; in particular round a path of infinitesimal size the 
integral also vanishes, i.e. curl Y s is zero* 

Because of the relation between V^ and S t Y s is called 
a lamellar or scalar potential vector field t S being its potential. 
Since the curl is zero, the field is also called a non-curl 
or irrotathnal field, 

2. The Inverse Square Law- Point Sources. In 
practice the greatest interest is found in lamellar holds 
based upon the law of inverse squares. Consider a point 
source of vector flux, around which the level surfaces are 
concentric spheres ; the flux lines are radial straight lines 
diverging from the source. At unit distance from the 
source the spherical surface has an area of 4Z ; if the 
source is of unit strength this is indicated by drawing one 
flux line per unit of area to represent the magnitude of the 
vector field at unit distance, A total vector flux of 47 1 
tines cuts normally through all spherical surfaces and hence 
the magnitude of the vector field at any point varies with 
the inverse square of the distance from the source, since 
the spherical area varies directly as the square of the 
distance. If the strength of the source is q y the flux is 

In some fields the point is an actual source of material flux, 
e.g* when liquid is entering at a point within a fluid ; the vector 
is then the outward radial velocity of the fluid. In most cases, 
however, the flux is non-material, e.g, of heat, or of gravitational 
or electric forces. Note that in the case of gravity the point mass 
can only attract another mass placed in its field; the force is 
essentially radially inward. But in hydromechanics, heat, elec** 
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tricity, and magnetism both positive and negative sources, the latter 
called sinks, can exist, making possible both outward and inward 
fluxes. Fluid or heat may either enter or leave; electric or mag¬ 
netic forces may be either repulsive or attractive. Masses act in 
the vectorial sense in the same way as sinks. 

It will be convenient to fix ideas by examining an actual 
example, namely, the field of electric point charges in 
vacuo. Two such charges of positive electricity have been 
shown by Coulomb’s experiments to be repelled with a 
force hq x q 2 /r 2 in the line joining them. If unit charge is 
defined as repelling a similar charge 1 cm. distant with a 


£ 



Fig. 28.—‘Electric Field of Point Charge 

force of t dyne, we have the ordinary electrostatic system 
of units, k - 1 and £j£*/r a is the force in dynes between 
charges and q 2 e.s.u. Let a positive charge q be put 
at O, Fig. 28, then the force per unit positive charge 
placed at P is called the electric force or field strength ; 
its magnitude is q/r 2 and its direction that of the unit 
radius vector r*, i.e, of r, then 
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The potential difference between two points A and B is 
defined as the work done against the forces of the field 
in moving the unit positive charge from A to B by any 
path. Since like charges repel* positive work must be 
done against the field forces to bring the charges nearer 
together. The work that would be done by the field forces 
is Edr ; hence the potential difference is, by definition, 

Now let A be at infinity—where £ is zero—and make B 
coincide with P ; then the integral becomes 



This is the potential at P, namely, the work done on unit 
positive charge in bringing it from infinity to the point P 
by any path. Using Equation 5.11 shows that 

£ - Jfi - - grad 

A similar argument can be adapted to other fields. In 
general, let F be a flux vector at a distance r from a source 
of strength q ; then the Newtonian potential is defined as 

4 > - q/r ;.( 7-0 

and F - - grad <f> -grad (q/r), . . (7.2) 

of which the three cartesian components are respectively 
F x - - dfifdx t F y - - dj>fdy and F t - ~ d^/dz . (7.3} 

The vector F points in the direction of the greatest rate 
of decrease of 0 , For example, if ^ is temperature it is 
clear that the flow of heat will be in the direction in which 
the temperature falls most rapidly. 

Since potentials are scalar functions of space they are 
arithmetically additive. Thus, if there arc several point 
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sources q lr q $ , . . at distances r u r % , . . from a point 
in the field, the potential is 

^ " [r| * + * * ’] ’ * * (74) 

This may be extended still further. If the sources 
are continuously distributed with density a per unit 
area over boundary surfaces the summation becomes the 
integral 



taken over all the surfaces. Each surface element is 
regarded as a point source yielding a Newtonian potential. 

In Equation 5.12 it has been shown that except at a 
point occupied by a source the field F is solenoidal, i.e. 
div F - 0 ; then 

. . . . . (7.6) 

he. Laplace’s equation is satisfied by the potential at all 
source-free points in space, stnee div grad $ is, by Equa¬ 
tion 5.12, zero, A field of this kind is, therefore, known 
as a Laplaaan field and <f> is called a harmonic function. 
Equations 7.1, 74 and 7.5 may be regarded as solutions 
of Equation 7.6 for various physical conditions. 

3 . Volume Distributions. In many problems the 
field is due to volume distributions of sources, each element 
of which acts according to the inverse square law. Fam iliar 
examples are the electric field due to a space charge of 
electrons or of ionized gas molecules, the attraction of 
solid bodies, the flow of internally generated heat (as in 
an electrically-heated conductor), &c. Let o be the density 
of sources per unit volume within a closed surface $ ; 
Q will be a scalar function of position. Then if dv is a 
volume dement in Fig. 29, @dv is an elemental source. 
Two cases arise and will be separately examined. 
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Let P be outside the surface s, then the total potential 
at P is 



taken over the volume enclosed by s. Also, since gdv is 
constant for each differentiating operation, 

F « - grad <j> and V 2 <£ s o, 

exactly as before. The field external to a volume distri¬ 
bution is solenoidai, he. the potential is Laplacian, 



Fig, at).—Potential of a Volume Distribution 
If P is inside the surface s there is a small difficulty, 
since r in the denominator of the integrand can now become 
zero. Surround P by a small spherical surface of radius 
e ; then, taking the volume integral through the space 
confined between s and the sphere, P is outside this region 
and the integral is finite. Now let the sphere be made 
vanishingly small; then in the limit, its volume approaches 
dv , which varies with Since r varies with £, dv/r 
depends on e 2 and vanishes as a becomes zero. Hence 
Equation 7,7 expresses also the potential inside the surface 
if it is regarded as the limiting value reached when a small 
sphere round P is reduced without limit * the held is not, 
however, Laplacian. In Equation 6 a it has been shown 
that wffien V •* grad (i/r) the total normal flux of V 
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through a surface enclosing a source of unit strength is 
- 4Jr independently of the position of the source. Hence 
when F - - grad <f> t where is due to sources of total 
amount q t no matter how distributed within the surface s t 
the flux is + 4 zrq. In the present case q is the volume 
integral of the source-density, so that 



n*F ds - 47E 



Applying the divergence theorem. Equation 6.1, 


J j" n-Fds “ j [ v 



so that div F - - - div grad $ 

and V 2 <f> - - 4^0, (7.8) 

which is Poisson's equation. Hence the field in a region 
where there is a volume distribution of sources is not 
solenoidal. There is divergence, be. total normal flux per 
unit volume, at any point, of an amount equal to 4jr times 
the density of sources per unit volume. 

It is obvious that the same results w r ill hold for two or 
more distributed sources, interpreting the integral. Equa¬ 
tion 7.7, as extending over all the enclosed volumes. In 
the space between them the field is solenoidal and the 
potential Laplacian * within each volume the field has a 
Poissonian potential and has divergence. 

4 . The Potential Operation. The integral, Equa¬ 
tion 7.7, taken over all space may be considered to be the 
solution of Poisson's differential equation, Equation 7.8, 
for any distribution of sources throughout space. The 
operation of finding the potential due to a given distribu¬ 
tion of source density g is indicated by the notation pot g t 
due to Gibbs. Then, making use of Equations 7.7 and 7.8, 

* " JJ jv - P ot e - - ~ pot W . (7-9) 
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Hence the operations 

“^P«( ) and V*( ) . . (7-9*) 

acting on a scalar operand are inverse, one undoing the 
effect of the other in much the same way that integration 
and differentiation are inverse operations* Applications 
of this notation have been much developed by Gibbs, 
and a few simple extensions of the idea will be found in 
the next Chapter* 

5 * Multivalued Potentials* In mathematical text¬ 
books it is shown that the operations grad, div, and curl, 


Y 



Fig* 30*— Multivalued Scalar Potential 


as well as theorems such as those of Gauss and Stokes, 
are valid provided that the scalar or vector point functions 
to which they are applied remain finite, free from discon¬ 
tinuities and single valued within any given region of 
space and that the same conditions also restrict their 
derivatives in any direction up to the second order (so that 
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operations such as V® may have a meaning). Fortunately, 
the conditions of finiteness and continuity are well observed 
in most fields of physics and engineering, but important 
cases do arise in electromagnetism and hydromechanics 
where the idea of scalar potential leads to multivalued 
functions* 

Consider first a simple mathematical example, the 
potential given by 

^ « arc tan (yfx ) - 0 ± for, 

where k is any positive integer including zero. Then the 
1 equipotential' surfaces are radial planes 6 - constant 
intersecting along the Z axis, as shown in Fig, 30(G). 
This intersection violates the condition of single-valuedness 
mentioned on p. 12 ; the axis is termed a singularity. 
Taking the gradient, 

grad # - V# - + j|; + “ 5) 



a X J - -(-sin 01 + cos 0 J) 

^ y g' 31 


The lines of the vector, grad are in the direction of the 
counter-clockwise tangent at any point P t perpendicular 
to the * equipotential * plane through that point* To cal¬ 
culate the line integral of grad $ we need 




xdy — ydx 
x 2 + y ' 1 


Consider a circular path round O in the plane of XY; 
then, as a is constant, from 

x - a cos 0 and y — a sin B 

we have 


dx = - a sin 0 d 0 t dy - a cos OdO and xdy - ydx m a-dd ; 


hence 
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along any portion of a circular path* Since it is inde¬ 
pendent of a the integral applies to any path,* For a 

closed path such as A t Jdfl is clearly zero, since to trace it 

out the angle subtended at O by the area is passed over 
once positively and once negatively by the radius vector 
from O. This is true even when is infinitesimal and 
hence curl grad <f> is zero, as required by Equation 5.3* 

Now trace a semicircular path ab, giving a line integral 
of 71. Likewise, a line integral along the semicircle a'h* 
gives - rr f i.e* the same magnitude but of opposite sign* 
But the beginnings and ends of the two paths are the 
same ; hcncc the value of the line integral depends on 
which path is followed as well as upon its end points. 
Round any closed path linked about the axis of Z the 
line integral is ± 2jz per tour, the upper sign relating to 
rotation in the sense from X to Y* Such a closed path 
cannot be reduced to zero without cutting through the 
Z axis and cannot satisfy Equation 5,3, The potential at 
P t defined as on p* 70, is multi-valued, being increased or 
diminished by 2 n for every tour made round the axis of Z 
before arriving at P from an infinite distance* 

Mathematicians describe a region of space in which 
these phenomena occur as * cyclic * about a singularity in 
the region ; in our problem the singularity is the axis of 
Z where all the level surfaces intersect. A path such as A 
which can be contracted to zero without enclosing or 
cutting through the singular point or points is * re¬ 
ducible 1 j a path such as abb* a* is 1 irreducible * since it 
cannot be made zero, but can only be made to fit more 

* Any path can be resolved into three components, viz* parallel 
to OZ t radially from OZ and around it. The first two can add 
nothing to the line integral, since the lines of flow of grad 4 are 
circular. At any radius a the path around the axis is add t for a 
•mall angular rotation. The magnitude of grad being i/n, 
the element of line integral is dO. 
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and more closely round the singularity* A cyclic region 
can be made acyclic by inserting impassable barriers to 
prevent the drawing of irreducible patlis, thereby isolating 
the singularity* For example, in Fig. 30(0) take the semi¬ 
infinite plane of ZX to the right of the Z axis as a barrier, 
shown shaded* Then 0 cannot now exceed zti and ^ is 
thereby made single valued and the conditions of potential 
theory are satisfied* 

The problem just examined Is easily seen to represent the 
magnetic field of a current f electromagnetic units in a thin straight 
wire extending to infinity along the Z axis if we take 

$ « — 21 arc tan iy/x). 

Then the magnetic force is 

II - — grad $ - (2t/fl)C— sin & 1 + cos & J) 

and the line integral of H round a closed path encircling the current 
is 4the so-called magnetomotive force. Remove the current, 
and let the half-plane of ZX to the right be occupied by a thin, 
uniformly-magnetized sheet of moment i per unit area with the 
north-polar surface uppermost (Fig, 30(6)), Then the magnetic 
field of the sheet and the current are identical in distribution, but 
that of the sheet is acyclic since no path can pass through it to 
encircle the wire* Its scalar potential is £ — zi(ir — 0 ), # being 
limited between o and zw (p. 90)* 

The fundamental axiom of potential theory is that the 
line integral of the vector field is independent of the path 
traversed between two points ; the theory is not applicable, 
therefore, without some artificial aid to problems which 
do not satisfy this condition. The most important case 
Is the essentially cyclic magnetic field of an electric current 
flowing in a linear circuit of any shape* Exploring the 
field, in the usual way, with a unit magnetic pole which is 
moved about from point to point, will result in changes in 
the amount of flux from the pole which links the circuit, 
inducing electromotive forces therein ; if the current is to 
remain constant, this will require adjustments to be made 
to the voltage of the battery in the circuit* For every 
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complete tour made by the pole round the current, potential 
energy ±4^*' must be credited to the pole; this energy 
is, of course, derived from the battery. Physically, there¬ 
fore, the scalar potential of the field must be multivalued, 
and its cyclic character has this physical explanation. To 
satisfy mathematical requirements, however, we replace 
the field of the current by that of a thin, normally-magnet¬ 
ized sheet with the circuit as its boundary edge—the 
equivalent magnetic shell. The potential due to the shell 
is m, to being the solid angle subtended by the circuit at 
any point in the field, and is single-valued as required by 
potential theory, p. 1 z; the field is then acyclic and the 
problem has been reduced to one in magnetostatics ; see 
further, p. 86. It is, however, not in agreement with the 
essential physical nature of a current field, and must be 
regarded as a mere mathematical convenience. Agreement 
can be secured by treating the field of a linear current as a 
degenerate case of the vector potential field of a current- 
carrying medium, as on p, 84. 


CHAPTER VIII 


THE VECTOR POTENTIAL FIELD 

1. The Magnetic Field of a Steady Current, 
Vector potential fields occur most commonly in electro¬ 
magnetism and in hydromechanics. It will be convenient 
to introduce the subject by working out the properties of 
the magnetic field of a steady current as a practical illustra¬ 
tion of the fundamental ideas; the reader will readily 
adapt them to any analogous physical problem in which 
he is interested. 

We have seen on p. 77 that the space round a circuit 
carrying a current is a cyclic region ; the physical meaning 
of this is that the work done on a unit magnetic pole in 
transporting it once round a constant current of i electro¬ 
magnetic units is ± 4JE1 ergs. This is the well-known 
circuital theorem . 

Consider now a medium, such as the interior of a copper 
conductor, which is carrying a current distributed in any 
given way. The current density at any point is c electro¬ 
magnetic units per sq. cm. of a surface perpendicular to 
the direction of flow. Draw any closed curve in the 
medium and take any open surface having the curve for 
its bounding edge. If n ds is an element of vector area on 
the surface the current through the element is the scalar 
product c*n ds, and that through the whole surface is the 
integral of this taken over the entire area. Now calculate 
the tangential line integral of the magnetic force II round 
the closed curve, in the sense of a right-hand screw relative 
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to the unit surface-normal n, as in Fig. 31; then the 
circuital theorem makes 


(jjH'ifl = 4nJ fc-n ds* 


Apply Stokes’s theorem, Equation 6.3, to the left side, then 

J jn*curl II ds “ 4JZ j* Jc*n ds t 

and for this to be true of any surface vve must have 

curl H - pic , (8.i) 


open surface s 



Fig. 31.—Curl of Magnetic Field due to Electric Current 

Physically this means that we put a very small plane area 
at a point in the medium and turn it about until the line 
integral of H round its boundary is the greatest possible; 
the amount of the integral per unit area is the magnitude 
of curl H and the positive normal points in the direction 
of the curl. But the equation relates this to the current 
density c, which must, therefore, be passing normally 
through the small area when it is in the orientation giving 
the greatest value of boundary line integral. 

It is to be noted that since div curl H is zero. Equation 
5,6, div c is zero also. Hence currents are solenoidal, 
he. flow in closed paths. Further, since there can be no 
such thing as a space distribution of free magnetic poles, 
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div II is zero also ; the lines of magnetic force must, 
therefore, form closed curves. 

2 * The Vector Potential* On p. 55 we have seen 
that when a vector has curl and no divergence it may be 
represented as the curl of another vector, called the vector 
potentials The reason for this name will now be explained* 
hirst take at any point in the current-carrying medium, 

II - curl A,.(8,2} 

where A is the vector potential of II. Using Equation 5.7 
on the above, 


curl H *= curl curl A - grad div A — S? 2 A - pic* 


We may, without any loss of generality, assume that A, 
like its related vectors II and c, is solenoidal; then div A 
is zero and 


V 2 A - - pic t . , . . {8.3} 

which is a vector form of Poisson’s equation relating A 
and c within the medium. Using Equation 54, put 

A • A x i + AJ + A z k and c - c z i + cj + c z k, 

then 

V 8 4 * - - 4 nc xt V 2 ^ = - ptc v and S/ 2 A s - - 

If r is the distance from an element of volume, where the 
current density is c, to any point in the medium at which 
the components of A are to be found, analogy with Equa¬ 
tions 7.8 and 7,5 gives 



and A t 



(8.4) 


Recombining the components, 


A 




pot V 2 A; * 


(8-5) 


the integrals extend throughout the volume in which c 
is distributed. 

Equations 84 show that the components of vector paten- 














8a 


VECTOR ANALYSIS 


dal are calculated from the components of c in exactly 
the same way as the scalar potential is found from a scalar 
source distribution of density g* Hence the name vector 
potential* Interpreting Equation 8*5, at the point where 
the vector potential is required draw a vector parallel to 
the current density c at any other point in the medium, 
the length of the vector being proportional to the magni¬ 
tude of c multiplied by dv and divided by r, the distance 
between the points* Repeat for every element of volume 
into which the region containing c can be divided and 
sum the vectors obtained ; i*c. integrate their effect at the 
point for the whole volume* The resultant vector is the 
vector potential at the point; the same result is more 
practically obtained by calculating A xt A yt A t and vectorially 
adding these rectangular components* 

By exact analogy with the scalar potential field, p* 72, 
the vector potential at a point outside the medium con¬ 
taining c is found by appropriate evaluation of the integrals 
Equations 84 and 8,5, but now 

curl H - curl curl A - V 2 A s o 

at the point Each component of A satisfies Laplace's 
equation* II is now a lamellar vector and can be expressed 
as the negative gradient of a multivalued scalar potential 
Q say* Hence at points outside a current-carrying medium 
the magnetic force is given by 

H - - grad Q - curl A* 

At the surface of separation which divides the current- 
carrying medium from outside space, the two expressions 
for A, in the curl and non-curl regions respectively, reduce 
to equality. Their analytical forms in the two media are, 
however, quite different** 

# A variety of practical cases will be found in the author** 
treatise, Electromagnetic Problems in Electrical Engineering (Oxford 
University Press, 1939). 
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3* The Potential Operation* The quantities A, H, c 
are all solenoidal fluxes and are related one to the other by 
the process of curling, Equations 8*1 and 8*2, and the 
calculation of potential, Equation 8*5* Summarizing these 
connexions, 

A - pot c 

H - curl A « curl pot c 

C « — curl H « — curl curl A — —curl curl pot C* 
4^ 4 JE r 

The fortuitous 1/4^1 is a consequence of the system of 
units that has been used ; it may be abolished by changing 
to the Heaviside rationalized system* The third relation 
shows that the operation denoted by curl cnrl/^z exactly 
annuls the operation of computing the vector potential. 
Again, the first and third shows that 

A - pot c - — pot curl curl A, 

so that pot/4Ji and curl curl also annul* Hence it follows 
that pot and curl curl are commutative operations ; we may, 
therefore, find the vector potential either before or after 
the operation of double curling* 

Although the three vectors considered are the only 
ones with physical meaning, we can extend the process in 
either direction* Adding a fourth, D - curl c, it follows 
from the construction of the table that H - pot D* Any 
vector in the continuous series is the curl of the one imme¬ 
diately above it (neglecting the factor) and the potential 
of the vector two below it. Now 

H « curl pot c - pot D - pot curl c j , (8.6) 

hence it is immaterial whether we calculate the vector 
potential before or after calculating the curl. 

These relations and many others have been developed 
by Gibbs into a complete calculus for manipulating sole- 
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noidal vectors, but to pursue it further would be beyond 
our present object* 

4 * Linear Currents* In electrical theory there is 
considerable practical interest in the magnetic field of 
currents flowing in infinitesimally thin conductors forming 
dosed circuits of any configuration* Such circuits are 
called linear, since the current is caused to flow along 
straight or curved lines, and are a close approximation to 
practical circuits of thin wire* In hydromechanics, such a 
linear curve constitutes a vortex filament in an incom¬ 
pressible fluid, the fluid being in rotation about the line 
as axis* In both cases the field is cyclic but non-curl* 
The electromagnetic case is often treated magnetostatic- 
ally by the * magnetic shell * conception, but we shall now 
see how vector analysis enables us to express the field in 
terms of vector potential for a medium free from current 
everywhere, except in the filament constituting the linear 
circuit; precisely similar arguments may be applied to 
vortex filaments in a fluid or to any other linear singularity* 

Let i be the current in the circuit and a its vanishingly 
small cross-sectional area* If t be unit tangent at any 
point of the circuit, we may put c - it/a and dv - a dl, 
where dl is an element of length of the circuit as shown in 
Fig* 32 \ then, 

c dv = itdl - i.dl 


and A" jlr .(8-7) 

will give the vector potential at a point P when the integral 
is evaluated round the closed circuit* The magnetic force 
at P is 


H - curl A »i curl 


<j)^ - f curl o l -dl « i(j) 


curl - dl, 
r 1 


since by Equation 8.6 the operations curl and pot are 
commutative* Using Equation 4*14, substitute i/r and 


THE VECTOR POTENTIAL FIELD 


*5 


t therein ; remembering that the curl is calculated at P 
and that t is not a function of that point, so that curl t 
thereat is zero, gives 

curl (t/r) - [grad (i/r)] x t - - (r t x t)/r 2 
from Equation 5.11, r x being unit vector drawn from the 
circuit element towards P. Substituting gives 




u sin 0 dl 


(S*S) 


r J r* 

where u is a unit vector in the sense of the vector product 



Fig* $z.—A mpere's Rule for Magnetic Force of Circuit Element 
(Biot-Savart Law) 

of Ti and t, Le. r lt t and u form a right-handed system* 
For example, if r x and t are in the plane of the paper u 
is perpendicularly outward toward the reader* Equation 
8*8 is known cither as Amp&re’s rule or as the Biot and 
Savart law * which states that the magnetic force due to an 
element dl of a closed circuit has a magnitude i sin 0 dl/r a 
and is normal to the plane containing dl and r ; the direc¬ 
tions of magnetic force, dl and r , form a left-handed system, 
since H is in the sense of — is. 

Calculating curl H, since Equation 5,3 shows curl grad 

* The general rule is due to Ampfcre for any shape of circuit- 
Biot and Savart f a proof applies only to straight conductors* 
















§6 VECTOR ANALYSIS 

(i fr) to be zero, the field of linear currents has no curl. 
As shown on p, 78, scalar potential theory can be applied 
to the cyclic non-curl field by replacing the current by a 
uniformly-magnetized shell having the circuit for its 
bounding edge. The north polar face of the shell is directed 
toward a point P, as in Fig* 33, when the current as viewed 


tfP 



Fio* 33.—Equivalence of Magnetic Shell and Linear Current 

from P appears to run counter-clockwise. If <o is the solid 
angle subtended by the shell, Le. by the circuit, at P t then 
the scalar potential is taken as 

<f> - m> t giving II « - grad <f> « — 1 grad 01 , (8*9) 

We now proceed to show that this agrees with the vector 
potential. 

Give the point P a small vector displacement dp in a 
sense to increase r, the vector joining P to an clement <fl 
of the circuit ; this is equivalent to fixing P and moving the 
circuit bodily in the opposite sense to the dotted position. 
The element sweeps out a vector area d\ x dp which 
subtends a solid angle Ti-dl x dp/r 2 at P, where r/ is 
unit vector in the sense from P to the element. The 
integral of this taken round the whole circuit is the total 
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change in solid angle subtended by the circuit when it is 
slightly displaced. Applying the rules for scalar triple 
products, Equation 2,21, gives 

~r( dl * rfp) - dp (^r * dj - pjp fjt x Jij, 

if Pi is unit vector along dp . Dividing hy dp and in¬ 
tegrating round the circuit gives 

dot f r/ w 
dp “ pl ’y^r * &: 

since P moves away from the circuit this is the rate of 
decrease of solid angle subtended by the circuit. But 
grad to is the total rate of increase of tu ; hence - p,-grad co 
is the component rate of decrease in any direction denoted 
hy Pi* equivalent to dto/dp above. Also if rj is unit 
vector drawn towards P from dl, r,' - - r t . Taking the 
scalar product with p, in Equation 8.9, 

HPi - - ip 1-grad a> = - - Pi-»<j>^^-^ dl. 

I'or this to be true of any displacement, 

H 

which is Equation 8.8, Hence the scalar magnetic poten¬ 
tial defined by Equation 8,9 gives the same magnetic field 
as the vector potential of Equation 8*7 : we can write, 
therefore, 

H - curl A « ~ grad <f>, 

where A « and <f >« proving the equivalence of 

the current and the shell in the field produced at any point, 
5 . Simple Examples of Vector Potential, The idea 
of vector potential is more readily grasped by considering 
simple examples of linear currents. First take the case 

7 
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of an infinitely long, thin wire carrying current along the 
positive direction of Z as shown in Fig 34(a). The lines 
of magnetic force are circles with their centres on the Z 
axis, the circles lying in planes parallel to the plane of XY ; 
in other words, since the wire is infinitely long, the field 
distribution in the plane of XY is precisely similar to that 



Fig. 34,—Vector Potential of Straight Current 


in any other parallel plane. Since all the elements of 
current are in the same direction, the vector potential at any 
point Pisa vector parallel to OZand can be found by a simple 
scalar integration. Writing dl - k dz and r 2 - z 2 + a\ 
Equation 8,7 gives 


A 



dz 

(c 2 + a*)* 



arc sinh 



« 2ik arc sinh(^/fl)J 

» 21k ’ log[# + \/(z 2 + a a )]| “ {C m - 21 log a) k, 


THE VECTOR POTENTIAL FIELD 89 

where C 9 is an infinitely great constant arising from the 
infinite extent of the 1 circuit * round which the integral 
has been taken. Since we are concerned with derivatives 

r 



Fig. 35,—Vector Potential of Circular Current 


of A, the presence of this constant is immaterial ; without 
loss of generality we can take 

A - A z k -1 - 21k log a m - ik log (jt* + 

Since A x and A^ are zero, the usual expressions for the 
components of curl A give, from Equation 4.9, 

//* - 'dAJdy - - ziyf{x % + y 2 ) *• - (si/a) sin 0, 

II y - - dA z /dx - 2 ix/(x* + y 2 ) - (21/a) cos 0, 

H t - o. (cf. p 75.) 

Surfaces of constant vector potential arc coaxial cylinders ; 
w'ith equal increments in passing from one to the next 
their radii increase in geometric progression. Lines of 
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vector potential are, therefore, closely packed together near 
the wire and fall in density at greater distances as shown 
by Fig- 3 #)- 

The corresponding multivalued scalar potential is easilv found* 
As in Fig. 30 (A) take the right-hand half plane of ZX as the 
equivalent shell (p. 77). Draw a sphere of unit radius with centre 
at P ; then the half plane subtends an area 477(7* — 0)/2*? — 2(n - 8 ) 
at the spherical surface. This is the solid angle subtended by 
the plane and Equation 8.9 gives 4 “ 21(77 — 8 ) as on p. 77* 
The components of — grad 4 * he. of II, are identical with those 
of curl A on p. 89, as the reader should verify by using Equa¬ 
tion 4.3 and 8 « arc tan fy/tf). 

Another simple case is the current in a circle of wire, as 
shown in Fig, 35, the circle lying in the plane of YZ. 
At any point on the axis of X the vector potential is zero, 
the point being equidistant from all elements of the 
circle* For a point off the axis, nearer therefore to one side 
of the circle, it is geometrically obvious that the resultant 
vector potential due to the whole circle is a vector parallel 
to the tangent at the nearest part. Let the point describe a 
circle about the axis, parallel to the current j then by 
symmetry this argument applies to any point on this locus* 
Hence the lines of vector potential for a circular current 
filament are parallel circles with their centres on the 
common axis. To give an expression for A and hence for 
H requires the use of elliptic functions and is beyond our 
present discussion. 

The diagram, Fig. 35, represents also in one meridianal 
plane the lines of flow of fluid about a circular vortex 
filament; the complete picture is obtained by revolving 
the diagram about OX * 


CHAPTER IX 


THE ELECTROMAGNETIC FIELD EQUATIONS OF 
MAXWELL 

I * General* The most striking example of the success 
of vector methods is provided by the theory of the electro¬ 
magnetic field* Although Maxwell introduced the funda¬ 
mental vector ideas into the theory, he did not make any 
use of vector analysis in manipulating his equations, relying 
entirely upon the usual lengthy cartesian methods. It is to 
Heaviside, Lorentz and many later workers that we must 
turn for the exclusive adoption of the concise vector 
equations and their treatment by purely vectorial methods. 
Our present object is to state the necessary equations 
and to use them as an exercise in the application of prin¬ 
ciples established in earlier chapters. We shall not, 
therefore, be concerned with questions of a purely physical 
nature, for the discussion of which the reader is referred 
to the special works mentioned in the Bibliography * 
One main assumption will be made throughout this 
chapter, namely, that the medium sustaining the field is 
homogeneous, Le. composed of one kind of material in 
the region considered, and also isotropic, i.e, having 
identical physical properties in all directions at every point, 
2 . Maxwell's Equations, At any instant let c be 
the current density and II the magnetic force at any point 
in the field, both being expressed in electromagnetic units. 

* See particularly F. W, G, White, Electromagnetic Waves 
(Methuen, 1934,) in this series of Monographs* Also see p. iiS, 
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Then by the method of p. So the circuital or magneto¬ 
motive force theorem, Equation 8.1 gives 

curl H - 4^eC. 

In deducing this relation, however, a purely conducting 
medium was assumed, c being the actual flow of electric 
charges per second normally through unit area. Maxwell 
obsen'ed, however, that electromagnetic phenomena take 
place in dielectric non-conductors, even in vacuous space, 
in which electric flow is impossible. To account for this 
he postulated a vector D, known as the electric displacemetit 
at any point and having the physical dimensions of charge 
per unit area, the time rate of change of which is called the 
displacement current density at the point. He regards the 
electric force £, which is the force per unit charge tending 
to move electricity at any point, as related to the dis¬ 
placement D In much the same way as stress is related to 
strain in an elastic solid. Using electromagnetic units, if 
k is the dielectric constant of the medium, 


D“#c£/ 4^.(g*i) 


In general a semi-conductor will exhibit both phenomena, 
i.e. true flow of electricity, the conduction current, together 
with displacement or, as an engineer would say, charging 
current. Writing the total current density as 


^ t?D k d£ 

€ « c + - 7 - c + — — t * 
at Apt at 

the first equation must be amended to read 


{9.2) 


curl H - 4J1C - 4jt(c + £)) -* 4J1C + «£, 


using the dot notation as an abbreviation for time differen¬ 
tiation, Maxwell’s assumption that a displacement current 
gives rise to a magnetic force in the same way as a con¬ 
duction current was a stroke of genius that led him to the 
theory of electromagnetic waves in space and the electro¬ 
magnetic nature of light; his conclusions were not veri- 
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fled by experiment until many years later. The total 
current flows in a closed circuit, i.e, it is solcnoidal and 

divC< *°. (9 3 ) 

Again draw in the field any closed curve and calculate 
the line integral of the electric force round it. This is the 
electromotive force round the curve and by Stokes’s theorem 
can be transformed into the surface integral of curl £ taken 
over any open surface with the curve as its edge. But by 
Faraday’s law the electromotive force is the rate at which 
the magnetic flux through the surface is decreasing. Let 
B be the magnetic induction or flux density, ft 

being the permeability, then 


<£ E dl - f 

f ri'curlC ds - - f- f f 

J J 

J *JJ 


which to be true for any surface bounded by the curv<> 
requires 

curl £ - - 6 - - ftil. 

The lines of magnetic induction are necessarily solcnoidal, 
i,e, 

div B ** pt div H — o. 

Summarizing, Maxwell’s equations for the field are 

curl H = 4J1C - 4^rc + *£, , . * (9.4) 

curl e - - ft - - yi ft ; . . . ( 9 , S ) 

together with C - c + £>,. (9.2 bis) 

div B - o, 
and div D » g T 

9 being the volume density of electric charge at any point. 
These general equations will be returned to in Section 5, 
3 , Conducting Media, When the medium is a con¬ 
ductor there is no displacement current and the main 
equations become 

curl H = 4 etc and curl £ -• - 6, 
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If y is the electrical conductivity of the medium* then 
c = y£, by Ohm's law ; the equations can now be written 
curl B ~ pipe and curl c - — yfi. 

Taking the curl of the first, Equation 5*7 gives 

curl curl B « grad div B — V - pift curl c - - 47^6, 

But div B is zero, so that the equation for B is 

* ... (9.6) 

Similarly, by curling the second equation and noting that 
div c vanishes, 

V 2 c “ pipyt .(9.7) 

These equations have been solved for many practical cases, 
e,g* eddy current induction in wires, bars or plates. 

They can be conveniently put in terms of vector poten¬ 
tial by defining 

curl A - B f * , * , * {9*8) 

together with curl B ~ piptc and curl c = - yfi* 

Now take the curl of curl A and choose A to make div A 
zero; then 

y*A - - 4^1/ic * • . • . (9.9) 

The third equation makes 

C - - yA.{9.10) 

Eliminating c 

V 2 A - ptuyh . . , . , (9.11) 

4 . Dielectric Media. In a dielectric there can be no 
conduction current, so that Equations 94 and 9,5 become 
simply 

curl H - and curl £ - — jw II, 

Curling the first, 

curl curl H - grad div H - V "H = k curl t - 

- Kti{dw/dt% 

Since div H - o, this can be written 

3 s 
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with precisely similar forms for B, £ and D, The 
bracketed operator is known as the dAlembertian, after 
the French mathematician d'Alembert, who first used a 
particular form of it in the theory of wave motion along 
elastic stringy. It is not difficult to show that Equation 
9,12 indicates that the components of II (and those of the 
other three quantities also) are propagated in the medium 
with a finite velocity 1/%/(*/*)■ I n ^ ree s P ace this is the 
velocity of light, a fact which led Maxwell to the conclusion 
that light is itself an electromagnetic phenomenon. 

5 . Energy Considerations. Calculation of the energy 
conditions in an electromagnetic field provides an excellent 
example of vector methods. Returning to the general 
case of a semi-conducting dielectric, consider an element of 
volume dv in which the electric and magnetic forces are 
£ and H respectively and the conduction current density is 
c. Then the Joule heat generated per second within the 
element is 

dP ~tcdv = (EjC* + EyCy + Ex c ) dx dy dz % 

the cartesian form enabling the reader to verify the vectorial 
statement. Now Equation 9.4 makes 

C “ — curl H “ £), 

4 n 

and dP - j”^- E-curl H - C t) J dv. 

Equation 4.12 gives 

div (£ x H) - H-curl £ - £-curl H, 

so that dP — 1 F—H-curl £ —— div f £ x H) — 

l_47i pi j 

Using the second main equation, 9.5, and also Equation 9.1* 
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But it is easy to see that Equation 3.3a makes 

H ' ft “ and £ ' £ " ( £i ). 

Also Gauss’s divergence theorem, Equation 6.1, trans- 
forms the volume integral 


JTf div * h > dv inti> J{ n ‘( £ * n ) *. 


the surface integral of the normal component of £ x H 
taken over the surface enclosing a given volume of the 
field* Finally 


* (913) 


which is Poyntings theorem. 

On the left side we have the sum of the Joule heat 
generated and the magnetic and electric energies stored 
per second in the volume. On the right is the total 
inward normal flux of a vector 


R “= (£ x H)/4^r t 
which is known as the Poynting radiant vector. It must 
represent the rate of flow of energy per unit area normal 
to R, Note that R is normal to the plane containing 
£ and H and forms a right-handed system with these 
vectors. Equation 9.13 states that the sum of the energy 
per second wasted in heat and that stored electromag- 
neticaUy in any element of volume of an e.m. field has 
entered that element as a flux of energy per second normally 
through its bounding surface. 


CHAPTER X 


ELEMENTARY PROPERTIES OF THE LINEAR 
VECTOR FUNCTION 

1. The Linear Vector Function. The simplest 
example of a linear vector function is provided by a vector 
U which is a constant scalar multiple of a second vector V, 
represented by 

V - AV, * * . . * (iQ.i) 

k being a constant. This simply means that U is k times 
the magnitude of V and in the same direction. If V is a 
point function of space, so also is U but altered in scale k 
times. Familiar examples are the relations between B 
and H or between D and £ in homogeneous, isotropic 
media. 

Rather more generally, let each component of V be 
multiplied linearly by a constant, different for each com¬ 
ponent. Using rectangular co-ordinates, if 

v - V x l + VJ + . . . (10,2) 

^d U- £ 41 + UJ+ UX . . * (10.3) 

then this more general type of linear vector function 
requires 

V x -= k r V XJ U v - k tJ V v and U z - k t V z . (104) 
The vector U is not only different in magnitude but also 
different in direction from the vector V. As examples 
may be mentioned the relation between D and £ in crystals, 
or that between stress and strain in a crystalline solid ; 
indeed any case of anisotropic media in which the rect¬ 
angular axes are principal axes of strain. 

Most general of all is the case of complete anisotropy, 
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where each component of U is a linear function of all three 
components of V. In rectangular co-ordinates the relation 
between the components can be written 

u * - + Ky v u + 

“ kyxVx + ky^Vy + f * . (tO.5) 

u* - k^v z + k^v v + k Q V t .) 

The nine coefficients characterize the transformation of 
the components of one vector into those of the other ; in 
general, transforms V$ into one of the three parts of 
* 4 * giving a rule whereby to remember the significance of 
the subscripts* The final result of the transformation is 

U - 0 V t .(10*6) 

tvhieh may be regarded as a generalization of Equation 
10,1 ; 0 is an operator turning V into U by a linear modi¬ 
fication of components specified by Equation 10.5* It is 
rather graphically described as a Cartesian tensor, or simply 
as a tensor , since the relation of strain to stress in an aniso¬ 
tropic elastic solid obeying Hooke's law is of this kind. 

The operator 0 is also sometimes called an affinor, for the 
following reason* Let V have the particular value of r, the 
radius vector from the origin to a point in space* Then U is a 
vector to a second point related to the first by a general linear 
strain ; the points are said to be in 1 : 1 correspondence or in 
affinity. 

2, Simple Types of Tensors, The essential part of 
a tensor operation is the array of coefficients, such as 
which in the general case are nine in number ; this can be 
symbolically represented, writing the array in the form of 
a matrix # , thus 



• More correctly a square matrix. The general type of matrix 
has mn dements in m rows of n columns. It represents the array 
of coefficients in m linear equations with n variables and is of a 
much more general type than we have space to examine here* 


LINEAR VECTOR FUNCTION 


99 


Tensors and their matrices are manipulated by certain 
rules, which are chosen to give results consistent with the 
algebra of groups of linear equations, such as Equation 10*5. 
A matrix is a convenient symbol for a whole group and 
enables the set of equations to be handled as a single 
entity* It is beyond our purpose to deal in any detail with 
matrix algebra, but one or two simple rules will be 
stated. First, if V - A + B we define 

0 (A + B) - 0 A 4 - ; 

he. the tensor operation is regarded as following the same 
law of distribution as for multiplication in scalar algebra* 
Second, if 0 1 and 0 % are two different tensors, 

0 X \ + $ S V - 0 % V + + <& 3 )V ; 

be* tensors follow the commutative law for addition* This 
the reader may readily verify by setting down the sum of 
two sets of equations like Equation 10.5, when it will be 
seen that the sum of two matrices is a new matrix in which 
the elements are the sum of the elements of 0 l and 0 %. 
Third, a tensor is regarded as negative when every coeffi¬ 
cient in its matrix is reversed in sign. 

If the columns and rows of a matrix are interchanged, 
the resulting tensor is the conjugate tensor , for which the 
matrix is 



When two tensors are such that one is the conjugate of 
the other, inspection of Equations 10.7 and 10,S shows that 

&xy "* kxi = k+x and k tJZ ™ kjy * * (10*9) 

If a single tensor satisfies this condition it is called sym¬ 
metrical or self-conjugate. It has only six independent 
elements and may be written 

Our work is limited to the cartesian tensor and its square matrix, 
for which m » n ■* 3* 
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r^jr Afj, Aj.j-1 

— I ^yy ^ya ■ • • 

' — ferr Af,,* Ah—I 


tym 


foyy ^gi 
k h 

^]fl ** 22 ~ 


(10.10) 


Symmetrical tensors are of frequent occurrence in physical 
problems and their meaning will be adverted to in Section 3. 

In particular when the relation between U and V is 
referred to principal axes all coefficients are zero except 
the diagonal k m k ZiJ corresponding to Equation 104* 

When two tensors are such that one is minus the con¬ 
jugate of the other, the coefficients must be such that 

*»-Ayy-Aa- **y-V> *»““*b and k r~~K {»•**) 

If a single tensor satisfies this condition it is called anti* 
symmetrical, anti-self-conjugate or skew. It has only three 
coefficients and its matrix is 


ro hnt hjx t 

o k yt L . , 

ha ~ hyz 0 ^ 


(lO*I2) 


Since a skew tensor shares with a vector the property of 
having only three components, the operation of 0^ on a 
vector V is exactly equivalent to the vector product of 
two vectors, since the final result is itself a vector U. 

3 . The Symmetrical Tensor* The symmetrical 
tensor, Equation io,ro, has a simple geometric interpreta¬ 
tion that will now be explained* Let 1 , j, k in Fig. 36 be 
unit vectors along rectangular axes X t F, Z. With the 
same origin let a, b, c be unit vectors along a second set 
of rectangular axes A t B , C. The two systems are related 
by the cosines of the nine angles between all possible 
pairs of unit vectors, one in each system ; these nine 
cosines may be represented by the matrix, 

i j k 

i\ 


*r 4 4 4 -i 

> m 1 m 3 m s 
:Lf!j u g n 9 J 
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e.g. 4 is the cosine of the angle between i and a, Ac. Since 
both systems are rectangular, it is clear that only three of 
the cosines can be independent, Le. one system can be put 
into the position of the other by three independent angular 
movements. Project the i, j, k system into the directions 
of a, b, c; then it is geometrically obvious that 

a ** 41 + 43 + 4 k 
b - 

C ® n x i + n J + fl^k. 


Y 



Now from Equation 24, a 3 - fo 2 — c s « 1 and from 
Equation 2a a*b = b e - ca - o ; also i 3 - J 2 - k a - 1 
and !■] - j-k » k*i * o. Calculating these scalar products 
gives 

4 * + 4 2 + 4 2 - »*i 3 + m * % + m * ■ n \ + n % + »a* - if 
l$n A + 4 * k » + 4^3 “ “ 

nJi + nj t + fi 3 4 * o f 

which are the well-known six relations between the cosines ; 
showing that only three can be independently chosen* 

If the relation of U to V is given in one set of axes. 
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what is it with reference to the other ? Let the a, b, c 
system be selected in such a way that they are principal 
axes ; then, as defined by Equation 104, 

U m - k a Vn U b - k b V b and U 6 - k c V< 

is the relation between the components along A, B and C 
respectively. Now resolve V„ V yf V tf the components of 
V in the X t Y, Z system, into the other system, thus 

V a - hV z + l M V w + hV B ) U m - kv 9 
V b - m t V s + V v + m y V x > and U b - k b V b 
V c - n t V x + n t V v + t H V t j U e - k c V t 

Now resolve U at U bt U e upon X t giving 
U x - l x U a + m x ll b + - (Vi 2 + 

Likewise resolving upon Y and Z, 

U v - (VA + + Kn&tW* 

+ (Vi* + V«s a + Mi a )^ 

+ (AdVa + K m t m * + K n % n xlV» 

U g - (kJJt 4* k^n ± m M + 

+(W»+^i«»+Mj«3)^+(W+^» 2 +^3^l' 
Comparing these three relations with Equations 10.5 it 
will be found that the coefficients satisfy Equation 10*9 ; 
the tensor relating U to V in the system of axes X f Y, Z 
is, therefore, symmetrical like Equation 10.10, It follows 
that any symmetrical tensor corresponds with a transfor¬ 
mation from the principal axes to another rectangular 
system of reference. 

It is possible to give a simple graphical construction for 
a symmetrical linear function. Using Equation 10.10 the 
function is 

u = 0 tvm v. 

Now project V upon U, i.e. calculate by Equation 2.10 
V U - V,U, + VyUy + V,U„ 
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Using the symmetry of Equation 10.10 in Equation 10.5, 

v-u - y- 0 tvn y - k a v K ' 4 k m v* + k u v* + 

AK v V t V v + k vt V y V, 4- kctVfVt] = 5, 

say, where 5 is a scalar point-function. But 

U* - kfxVt + k^Vy + k a V, « \ ■?— 

and similarly U v - \ and U, - J 

But these are the components of £ grad S in terms of 
V K , V y and V z as co-ordinates, Equation 4.3, so that 

U - i/jl 4- E/jj 4- U,\l - £ grad S ; . (10.13) 



Fio. 37.—The Tensor Ellipsoid 

it follows, therefore, that U is a vector perpendicular to 
the surface S - constant in the direction of the outward 
normal. But S ■* constant is an equation of the second 
degree in the rectangular components of V ; regarding 
these as co-ordinates defining the extremity P of the 
vector V, the locus of that extremity is a quadric surface 
or conicoid. In particular, if S - 1 the surface is, under 
certain conditions, the tensor ellipsoid, and 

V*U - (resolute of V in direction of U)t/ - 1 {10.14) 

Hence this resolute in the direction of grad S is the recip¬ 
rocal of the magnitude U of the vector U. Fig. 37 shows 
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this geometric relation for a tensor ellipsoid set in the direct 
tion of its principal axes, 

4 , Resolution of a Tensor. Return now to the 
general tensor, Equation 10,7, and apply the distributive 
law to it and its conjugate ; then 


0 - 1(0 + 0 f ) + 1(0 - 0 C ) , . (10,15) 

If the reader will write down the matrices for 0 + 0 r and 
0 _ 0 c by adding corresponding coefficients, it is easy to 
see that the former is a symmetrical and the latter a skew 
tensor. Hence any tensor relation between two vectors 
is equivalent to the sum of a symmetrical relation and a 
skew relation. The first corresponds to a change of axes 
without deformation of the system of reference, as has been 
shown, and the second to a transformation of the angles 
between the axes themselves. 

5 , Repeated Tensor Operations. In vector theory 
we frequently need to express relations between more than 
two linearly connected vectors in an anisotropic medium. 
Thus in electromagnetism there is a tensor relation between 
H and R, and a second relation between R and the pondero- 
motive force on magnetized solids. In a general way, if 


• and 


U = 0 X V where 0 X - 


V - 0 jW where 0 2 - 




what is the meaning of 

U - 0 X 0 3 W ? 


The answer is easily found by writing out the three rela¬ 
tions between the components of U and those of V, Equa¬ 
tion 10,5, and the analogous set giving the connexion 
between the components of V and those of W. From the 
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second set of equations substitute the components of V 
in the first set, giving U in terms of W, The final matrix 
for the complete transformation can be verified by the 
reader to be 

C ikjxfixz+kxzC&y (kxxCxy + k *4 ~] 

{hyxt-zx + kyyCyx + kyzCzz) (hyjfixy+kyyCyy + kyjfity} fcyyGyt + kyzCzz) I 

(kscCsx + ktgCpx +ex) ( A«Ca;y4- ^zyC {i y + ku^ty) (&&£&+ ktyCyi + 

An inspection of this apparently complex matrix will show 
that its construction follows a very simple rule : The 
element in the rth horizontal row and cth vertical column 
of 0 x 0 t is obtained by multiplying each element in the 
rth row of 0 X by the corresponding dement in the cth 
column of 0 2 and adding the results. It is usual for¬ 
mally to define 0 x 0 t as the 1 product * of the two matrices, 
the rule for the 4 multiplication * being that just stated. 

It is obvious that 0 < 0 X will, in general, be quite different 
from 0 x 0 2 , hence the rule for matrix * multiplication *, 
as also for successive tensor operations, is not commutative. 
The algebra of tensors and their matrices, therefore, fails 
in the commutative law for multiplication exactly as we 
have earlier seen to be true of vectors, 

A matrix of particular interest is 

fi 0 °1 

' t> - o i o L 

L° ° *J 

which, as the reader can verify, has the important property 
that 

v 0 - 0 v - 0 . 

Hence v has the same effect in tensor algebra as unity in 
ordinary scalar algebra, since it leaves a vector unchanged, 
and is called the unit tensor or unit matrix. 

6. The Dyadic. Although tensors are widely used to 
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express linear vector relationships in an anisotropic field, 
there is another method by which this can be done. The 
dyadic, introduced by Gibbs as an essential part of his 
treatment of vector analysis, is preferred by some writers, 
and its use will now be briefly explained. 

Returning to Equation 10.3 and inserting the values of 
U z , V p and U z given by Equation 10,5, the expression for 
U can be written as 

U -fei+VJ+ft. x k) v s+ (k xu ukj +^k) v v + 

(A J; i +kJ+k u k)V t 

- lV x +mV v +nV, 

say. Using Equation 2.9, V x - 1 -V, V v - j*V and V, - k*V, 
giving 

U-li V+mj‘V+nk*V=(li+mj+nk)‘V=y V . . (10.17) 

The bracketed expression is known as a dyadic and consists 
of the sum of three terms, li, mj, nk, which are its dyads. 
Each dyad is made up of two vectors in succession, the 
first vector being called the antecedent of the dyad (such 
as 1, m, n) and the second vector its consequent (such as 
1 , 1 , k). It is to be noted that a dyad is not a simple product, 
but merely a symbolic juxtaposition of two vectors. In 
interpreting the operator V 7 * the vector nearest to V in any 
dyad is involved in a scalar product with V, this product 
then being a scalar multiplier for the remaining vector of 
the dyad. Comparing Equation 10.17 10.6, we see 

that the operation of a dyadic on a vector by means of this 
scalar product rule is equivalent to the operation with a 
cartesian tensor. 

In Equation 10,17, we interchange W and V, 

U'«Vll+V mJ+Vnk-V (H+mj+nk)-V. . (10,18) 

which is a linear vector function different from U and called 
its conjugate. In Equation 10,17 the dyadic appears as a 
prefactor and in Equation 10.18 as a postfactor. 

Starting again from 

V - V x \ + V v m + V t n 
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and commutating the scalar products for V xt V u > V t gives 
the identity 

U»V-U+Vdm+V'kn-V<U+lm+kn)-V*!F c = F*V(io.i9) 
The dyadic W c , in which the rdies of antecedent and con¬ 
sequent are interchanged, is the conjugate dyadic of W, 
Hence any dyadic used as a prefactor is equivalent to its 
conjugate used as a postfactor. If a dyadic is such that 
iF-V *■ V*!P it is self-conjugate ; if 5 P-V • — V*W it is 
anti-self-conjugate. These definitions are exactly parallel 
with those for the corresponding tensors on p, too. 

A dyadic can easily be expressed in terms of the nine 
possible dyads of the unit vectors i, j, k, and in this so-called 
nonion form the relation between dyadic, tensor and 
matrix are even more clearly displayed. On p. 106 substi¬ 
tute in U for V xt V vt V t in terms of the unit vectors, then 

W - (hji + VJ* + kjd) 

+ OWJ + *w4J + K-M) 

+ (A^lk + /tyjk + A«kk) . . (10.20) 

and also 

^[(kJUkJhkJk) . 

+(*a»i | +*i™ij+*xjk) + (AttiJ+AyJj+Afjg) 

+(*JcMJtf+* <r kfc)J -f (kjk+k Jk+A„kk)] (10.21) 

which arc exactly analogous to the matrices, Equations 
10.7 and 10.8. 

If a dyadic is such that when applied to a vector, either as 
a prefactor or as a postfactor, it gives the vector itself, it is 
called the unit dyadic or tdemfactor. For this to be the 
case all coefficients Aa# for which tx + ft must vanish and 
kxx “ ■* Kz “ l * Then, as is obvious, 

/ - ii + jj + kk makes / V - V / a V, . . <10.22} 

which is equivalent to the unit matrix, since this also 
corresponds with U - V, The idemfactor and unit 
matrix occupy the same place in dyadic algebra and tensor 
calculus as unity in ordinary scalar calculations. 
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It is outside the scope of this book to enter into further 
details of the dyadic algebra developed by Gibbs or that 
of matrices used in dealing with tensors. It will be 
sufficient to say that both form complete, parallel and 
equivalent systems, and have been applied to a wide 
variety of problems in pure and applied mathematics. 
They are of the non-com mutative class in regard to certain 
products which, as with vectors, are matters of definition. 

7 . Application of Linear Vector Functions. It is 
a matter for personal preference whether one uses the 
tensor or dyadic notation in practice, though present 
tendencies would appear to show that tensors and matrices 
are nowadays the more usual. The simple cartesian 
tensors we have considered are fundamental in all prob¬ 
lems of anisotropy, whether magnetic, electric, thermal or 
elastic, and also in certain purely geometric applications to 
kinematics and kinetics where transformations of rect¬ 
angular axes are involved, as for example in the treatment 
of finite rotations. 

More complex tensors, which we have no space to 
discuss, occur in the transformation of space-time co¬ 
ordinates from one observer to another in the theory of 
relativity. Also, in the general theory of electromagnetic 
machinery the properties of a number of electric circuits 
in relative motion are connected by sets of linear equations. 
The use of tensor and matrix notation is a natural means 
for investigating such problems as Gabriel Kron has 
recently shown with great success, but this very special 
application is entirely outside the purely vectorial scope 
of this volume and belongs properly to the general algebra 
of tensors. 


POLAR CO-ORDINATES 


The most important operations of vector analysis have been 
expressed in the text in terms of components referred to rect¬ 
angular or cartesian co-ordinates. Summarizing, they are 

„.d s- vs - H» + Ifj + If... 

dlvV-7-V.<4.6, 4 8) 


*■»-(¥-©MS-E)! 

• ■ 




(4.9, 4.10) 
► * (5-1) 


S is a scalar and V a vector point function. These usually suffice, 
but in some problems it is often more convenient to use some 
other system of co-ordinates specially adapted to the matter in 
hand ; two of these will now be examined. 

When there is symmetry about an axis, such as OZ, the 
cylindrical polar or columnar co-ordinates shown in Fig. 38 
are appropriate. The position of a point P is specified by the 
polar co-ordinates r t 8 and the axial distance z t these being 
related to the cartesian system by 

x — r cos 0 , y — r sin 8 and z sa 3 . 

It is possible to transform the above operations from one system 
to the other by purely mathematical means, but it is more con¬ 
vincing to work out the operations afresh for the element of 
volume dr, rd $, d:z shown to a larger scale in die lower part of 
the diagram, (r t 8, z) being the co-ordinates of its mid-point. 
As an example, if V w be the radial component of a vector V at 
P the flux through the clement in the r direction is 

( r, + i^jpfr)(r + idr)<28 ds - (^V r - - J dr)d8 dx, 

which to the third order of small quantities is 
d$ dr dz + Vrd$ dr dz. 
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Since the volume of the element is r dB dr dz, the radial com* 
ponent of div V is 


5F, Vr 
dr r 




The other components are found by a similar process. Like¬ 
wise, by the method of pp. 43-5 the components of curl V can 


r 



be calculated ; this is left to the reader as an exercise. Let 

a, b, k be unit vectors at P, a being in the positive direction 
of r, b in the sense of the positive tangent to the circle of radius 
r at P (Le, along the arc rdO) and k along QZ. Note that a, 

b, k form a mutually perpendicular, right-handed system ; a 
and b, though of unit magnitude, are functions of the position 
of P, Then, if V, $ Ve, V ■ are components of the vector point 
function V in the a, b and k directions respectively, it is easy 
to show that 


grad S •* 
divv» 
curl V « 


div grad 5 — 
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„ e 05 I hs. t bS. 

7S *5;“ + ;55 b + 5; ki 

w)*' 


Y 



When there is symmetry about a point, as in many geometrical 
and dynamical problems, the spherical polar co-ordinates shown 
in Fig, 39 are useful. The position of P is given by the radius 

8 * 
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r, the angle of co-latitude between OZ and r, and the angle of 
Longitude ^ between the 20 X plane and that which contains 
02 and r, The relation to the cartesian system is 

x - r sin $ cos <f> t y - r sin 0 sin ^ and z - r cos 0, 

Taking (r, 6 , 4 >) as the co-ordinates of the middle of the volume 
element dr , rdd f r sin let a, b ( c be a right-handed orthogonal 
set of unit vectors at P in the positive directions of these incre¬ 
ments ; note that a, b, c are functions of the position of F, 
If V r , V$ t F* are components of a vector point function V in 
these directions, it is easy to prove for the volume element con¬ 
sidered that 


I as I a£ 

grad o — V jS — £ a + - t-i b -I- — : — —c: 

ar r d$ r sin 0 b^ 1 


“ ,v -” v -nb[s»'"-!'-' 5 ] 




^rdVe sin &) + Vt > 

&0 V ' r sin 0 B $ * 


+ r[ain 6 I* “ & (rp '*>] b + ~ w] C > 


divgrad 5 - V'S - ^ (r^—) + —- - (sin $*-) 

r 1 dr\ dr j r 1 sin & 50\ S1 bB / 


i b*s 

' r> sin* 0 


PROPERTIES OF V AS A FORMAL VECTOR 

On p, 33 a definition of V in cartesian vector form has been 
given, but it must be remembered that V cannot be regarded 
as a true vector since b/bx* b/by, a fbz arc differential operators 
and not true scalars. The symbol V actually denotes an in¬ 
variant vector differentiator and, as such, requires an operand, 
whether scalar or vector, upon which it can appropriately act. 
Its expression in vectorial form, therefore, can only be taken as 
symbolic or formal, inasmuch as V behaves like a vector in formal 
associative and distributive differentiations such as 
V{U + IP) - VU + VIP 
or V(UIV)-(VU)W + t/(V IP) 

where 17 and IF are scalar point functions ; these and many other 
simple operations the reader will find in the text-books and may 
easily verify. 

It must not be assumed, however, that V in more complex 
operations necessarily obeys all the rules of vector algebra, since 
V is only formally and not actually a vector. Care must be 
taken to see that rigid application of those rules does not destroy 
the different!ative effect of the complete operation. Particular 
care is needed in using the commutative law; examples are 
found in the operations V'V and V x V on a vector point 
function V, which become meaningless when commutated (see 
PP- 42 and 45). We conclude, therefore, that V is essentially 
and primarily an operator ; secondarily, and then only as an 
analytical convenience, is it expressed in cartesian vector form. 
In arranging more complex operations the object should always 
be to choose such a use of the vector laws as will leave differentia- 
tive operations available to act on the operand. 

An important practical case is that of the operator defined by 
C*y where 

c - C x i + CJ + C t k 

is a constant vector, Le. one independent of x , y , x, Taking the 
scalar product gives 


b 


b 


b 
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which is a scalar differentiator applicable cither to a scalar or to * 
vector point function. Note that V*C would vanish and thus 
have no meaning. In particular let 

c « cj + + cX 

be a vector of unit length in a given fixed direction defined by 
c xt c yt c St its components along the axes, i.e* by the direction- 
cosines of c, Then the operation 

a a a 

C - V - e *& + c ^ + c ‘Q? 

Apply this operator to a scalar point function S, then 


as 


as 


as 


(c-V)S-c^ +« lgs . 


Now from Equation 4,3 


as as as 


gmdS-VS-^i + ^J + a-k, 


is the greatest rate of increase of S at a point in magnitude and 
direction ; and c(grad S) is the resolved or component magni¬ 
tude of the rate of increase in the direction of e, a scalar point 
function. Calculating this scalar product and comparing results 
shows that 

(c* V)S - C (grad S) « c*(y S). 


This is known as the directional derivative of S in the direction 
of c, expressing the magnitude of the rate of increase of S in this 
specified direction at any point. Similarly, for a vector point 
function 


where 


V 

(c-v)V 

(C'V)Fi 


V^+VJ + VJl, 

(c'VWJ + (c'v)^J + {C'VWJk, 


BV* 

: dx 


BV* 

*h 


oV* 

dz 


is the directional derivative of V in the direction of C; it is a 
vector point function. 

This theorem provides an alternative proof of a result given 
on p, 47, namely, to calculate V * to * r, where 
<0 “ + cyjj + <uj , k 

is a constant vector and 


r - xi + yj 4 ark 
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is a variable one. Expanding as a vector triple product by 
Equation 3.22, 

V x to x r - {y r)io - (V’to)r 

But 


/ a & t \ f>x by &* 

v ' r - (fe + ^ + H;)< xl+yi + * k > - 5 + &? + 5 " 3 - 

Now from above V*to has no operational meaning ; we interpret 
it, therefore, as ujV. Substituting to for C, 


and 


w-v - 

/ B _ & i \ 

(W V)r - +. + yj + at) 


** tijji *f at v j + cu £ k — CO. 

Finally, 

V k w f" 3(1)- to** ata, 
as found by direct expansion on p. 47, 
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NOTATION 


Thb notation used in this book is that of Gibbs, with the additional 
adoption from Lorentz of the term * grad * to represent the 
differentiation of a scalar by the operator V ■ On the Continent 
certain differences of notation occur and these also appear in some 
English books (often translations). These are (i) the use of 
AB instead of A*B, or the employment of round brackets (AB) 
to denote the scalar product, (ii) The use of [AB] for the vector 
product AxB. (iii) The substitution of * rot 1 for * curl \ (iv) 
The writing of A, B or of A; B for the dyad AB. These differences 
are easily remembered and vector equations can be read with 
equal facility either in the Gibbs or the Continental system after 
a little practice. 

Heaviside uses an unsym metrical notation that makes hia 
writings unnecessarily hard to follow ; he uses AB and FAB for 
the scalar and vector products respectively. 

In reading vector equations, expressions such as A«B and 
AxB are read * A dot B * and * A cross B * respectively. 
Similarly, VS is read 1 del S \ V*V is * del dot V ' and V x V 
ta 1 del cross V 

The notation for line, surface and volume integrals follows 
Gibbs, i.e, 

J( W Jf< )di and JJJ( )dv 

respectively, the number of integral signs being a useful reminder 
that dl has one, ds two and dv three dimensions of length irrespec¬ 
tive of any particular co-ordinate system. 

Historically the terms used in vector analysis are of great 
interest. Although the operator V was introduced by Hamilton 
he left its development to his disciple P, G. Tait, who adopted 
the name 1 nabla * for this symbol, following a suggestion of 
W. Robertson Smith, Maxwell suggested * atled but the 
name * del * due to Gibbs is preferred for shortness and euphony. 
To Maxwell are due the terms 1 slope f (for 1 grad *), * conver¬ 
gence 1 (for — * div *), * curl 1 and 1 concentration 1 (for V *). 
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NOTE ON MAXWELL’S EQUATIONS 




The deduction of Maxwell’s equations for the electro* 
magnetic field given on p* gz et seq . follows that given in 
Chap. IX, Vol. II of his Treatise on Electricity and 
Magnetism, all quantities being expressed in electromagnetic 
units. It is very common practice, however, to state c, 
H, pL and B in the electromagnetic system, but £, k and D 
in the electrostatic system* If is the velocity of light 
in free space the last three quantities are expressed in 
electromagnetic units by writing v £, k/v % and B/v* 
Equation 9.2 then becomes 

C«c 4 -If>^cH -— t \ 

v 

and tlie field Equations 9.4 and 9*5 arc 

curl II “ arrC axc -f — £ 

V 

curl 

V v 

together with 

div B =* O 
and div D =* q 

The reader will find it a useful exercise to work out any 
modifications introduced into Equations 9.6 to 9.13 by 
the use of this mixed system of units* 
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partial, 32 

del operator, 33, M3 
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for scalar product, 17 
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Fields, scalar, 12 

level surfaces in, 12 
vector, 12 

lines of flow in, 13 
Finite rotations, 7, 46 
Flux, line of, 12 
total, 20 
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VECTOR ANALYSIS 


Gauss''a theorem, 59 
applied to inverse square law, 
61 

Gradient of divergence, 52 
Gradient of scalar field, 33, 35 
as operator V, 35 
meaning of, 35, 38 
in rectangular co-ordinates, 
36 

in cyl indrical co-ordinates, 111 
in spherical co-ordinates, 112 

Helmholtz's theorem, 55 

Idemfactor, 107 
Inverse square law, 57, 63 
point sources, 56 
Gauss's theorem, 61 
applied to electric field, 69 
Irrotational motion, 48 

Lamellar vector, 37, 38, 54 
closed line integral of, 38 
Laplace's equation, 54, 71 
Laplace's operator, si 
as div grad, si 
in rectangular co-ordinates. 
Si 

in cylindrical co-ordinates, 111 
in spherical co-ordinates, in 
Level surface, 12, 34, 67 
Line integral of vector, 18 
meaning of, 19 
in lamellar field, 37, 38 
Line of flow, 12 
Line, vector equation of, 7 
Linear currents, theory of, 84 
equivalent magnetic shell, 86 
Linear vector function, 97 
principal axes, 97 
cartesian tensor of, 98 
conjugate of, 99, 106 
skew, 100 
symmetrical, 100 


Linear vector function —contd 
resolution of, 104 
repeated, 104 
dyadic of, 107 
nonion form of, 107 

Magnetic field, of straight wire, 
77 

of current-carrying medium, 
79 

of linear currents, 84 
Magnetic shell, 78 

equivalent to linear current, 
86, 90 
Matrix, 98 
conjugate, 99 

skew and symmetrical, 100 
resolution of, 104 
product, 105 
unit, 105 

Maxwell's equations for e.m« 
field, 9i ( 118 
in conductors, 93 
in dielectrics, 94 
energy relations, 95 
Multivalued potentials, 75, 82 
cyclic and acyclic regions, 76 

Nabla operator, 33 
Newtonian potential, 58 
Normal, positive sign of, 19 
Notation, comparison of, 117 

Operator, del, 33, 113 
grad, 35, 36 
div, 40, 42 
curl, 45 

div grad, 51, 52 
grad div, 52 
div curl, 53 
curl curl, 53 
potential, 73, 83 
tPAlembertian, gj 
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Partial differentiation of vec¬ 
tors, 32 

Plane, vector equation of, 17 
Point function, 12 
Point source, 56, 68 
Poisson's equation, 54, 73 
Polar co-ordinates, 109 
Potential operator, 73 
Potential, scalar, 38, 67 
of point source, 57, 68 
Newtonian, 58 
in electric field, 70 
of surface distribution, 71 
of volume distributions, 72 
multivalued, 74 
Potential, vector, 55, 81 
in lamellar region, 82 
of linear currents, 84, 90 
Poynting’s theorem, 96 
Products of more than three 
vectors, 27 

Products of three vectors, 25 
Products of two vectors, scalar, 
14 

vector, 20 
Pseudoscalars, 26 

Rotation, 43, 48 

Scalar algebra, 2 
Scalar, definition of, 1 
Scalar potential fields, 12, 38 
of point source, 57, 63 
level surfaces in, 67 
lamellar nature of, 68 
of surface distributions, 71 
of volume distributions, 72 
multivalued, 74 
Scalar product of two vectors, 
14 

commutative property of, 15 
conditions for zero, 15 
■elf, 13 

application to unit vectors, 16 
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Scalar product of two vector* 
—eon id 

distributive law, 36 
in rectangular co-ordinates, 17 
differentiation of, 31 
Scalar triple product of vectors, 

25 

in rectangular co-ordinates, 

26 

Solenoidal vector, 41, 54 
Solid angle, 61 
and magnetic shell, 78 
Spherical polar co-ordinates,! 11 
Stokes's theorem, 63 
Subtraction of vectors, 6 
Surface integral of vector, 19 
meaning of, 20 

Tensors, cartesian, 98 
conjugate, 99 
skew, too 
symmetrical, 101 
representation by ellipsoid, 
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resolution of, 104 
product of, 104 
unit, 105 

Triple product, scalar, 25 
vector, 27 

differentiation of, 31 

Unit dyadic, 107 
Unit tensor, 105 
Unit vectors, 10, 16, 22 

Vector area, 22 

of closed surface, 23 
Vector, definition of, 2 

graphical representation of, 4 
scalar multiplication, 5 
polar, 5 
axial, 6 

Vector fields, 12 
classification of, 54 
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Vector potential, $$* 79 » ® l 
calculation of, 82 
in lamellar region, 82 
of linear currents, 84 
of straight wire, 83 
of circle, 90 

Vector product of two vectors, 

30 

noo-commutative property, 

31 

distributive law, 34 
in rectangular co-ordinates, 

*4 . , 

differentiation of, 31 
divergence of, 49 


Vector quantities, 1 
invariance of, 10 
Vector triple product of vectors, 
2 7 

Vectors, addition and subtrac¬ 
tion of, 6 

associative law for sum, 8 
components of, 8 
aum in rectangular co-ordin¬ 
ates, 11 

products of two, 14 
scalar, 14 
vector, 20 
triple products, 27 
diff erentiation of, 39 
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